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ABSTRACT 


Title of thesis: Meter-wavelength observations of pulsars 

using very-long-baseline interferometry 

Nancy Ruth Vandenberg, Doctor of Philosophy t 1974 

Thesis directed by: Professor William C. Erickson 

This thesis presents the results of an investigation of the angular 
structure imposed on pulsar radiation due to scattering in the inter- 
stellar medium. The technique of very-long-baseline interferometry 
(VLBI) was used to obtain the necessary high angular resolution (~0.1 
arcsec). The interferometers formed by the Arecibo (305 m) , NRAO (92 m), 
and Sugar Grove (46 m) telescopes were used at radio frequencies of 196, 

111, and 74 MHz during seven separate observing sessions between 1971 
November and 1973 February. The observations were concentrated on the 
Crab nebula pulsar (PSP. 0531+21), and several strong, low-dispersion 
pulsars were also observed. 

A crude visibility function for the Crab nebular pulsar was obtained 
during a 144-MHz VLBI experiment in 1971 March. The visibility function 
is consistent with a Gaussian brightness distribution having a half-power 
diameter of 0.0540.03 arcsec. Since only a small range of points in the 
(u,v) plane was sampled in the later experiments, the fringe visibilities 
were also interpreted in terms of a Gaussian model, yielding half-power 
diameters of 0.0340.01, 0.0740.01, and 0.18+0.01 arcsec at 196, 111, and 
74 MHz, respectively. A Gaussian distribution supports the view that the 
apparent angular sizes are caused by multipath scattering in the Interstellar 
medium. The apparent sizes obtained at all four frequencies scale as 



wavelength to the power 2. CKO. 4, supporting theoretical scattering laws which 
predict a wavelength-squared scaling. At 74 MHz the measured size agrees 
with that derived from an independent interplanetary scintillation experi- 
ment. The angular sizes were not observed to vary during the entire ob- 
serving period. 

The technique of differential fringe phase was used to show that 
the pulsar and the "compact source" in the Crab nebula are conincident to 
within 0.001 arcsec which corresponds to ~2 a.u. at the distance to the 
nebula. The total time-averaged flux from the compact source was also 
obtained. The fluxes measured at 74 MHz (105 Jy) and 111 MHz (40 Jy) lie 
on the extrapolation from higher radio frequencies of the spectrum of the 
pulsing flux, and the 196-MHz flux (8 Jy) agrees with the pulsar's flux. 

Variations in the total flux of up to 25% were observed. 

The data were also analyzed to yield the Fourier components of the 
Crab nebula pulsar's pulse shape by coherently de-dispersing the signal. 

To obtain the correlated pulse profile, the components (including the n=0 
harmonic) were added in a Fourier sum. At 196 MHz the reconstructed pulse 
shape is identical with the single-telescope profiles observed at Arecibo. 

The 111-MHz pulse profile contains both pulsing and non-pulsing components, 
each of which contributes approximately equally to the total flux. The 
reconstructed pulse shape agrees with the single-telescope pulse profile 
observed at Arecibo. Mo pulsing co’mponent was detected at 74 MHz. 

The ratio of pulsing to total flux (Ip/I^,) at 111 MHz was found to 

vary over the one and one-half years of observing. During the winter of 

1971-2 there was more pulsing than steady power (Ip/I T ~ 0.6) whereas in 
the winter 1972-3 the steady power was the dominant component (Ip/I^ ~ 0.4) . , 



The fringe visibility of the time-averaged pulsing flux and of each 
Fourier component is the same as that of the total flux. Both the pulse 
shapes and ratios of L^/I^ are similar for the unresolved and partially 
resolved fluxes. These statements show that there is no dependence of 
fringe visibility on pulse phase. This fact, and the observed combin- 
ation of temporal and angular broadening, rule out scattering theories 
which employ a simple single- screen model of the interstellar medium. 

A two-sc;:een model in which one screen is in the Crab nebula improves 
the agreement, but the model is still not consistent with all aspects 
cf the observations. A new model is proposed in which the temporal 
broadening is dominated by scattering in the small filament which lies 
in the nebula in front of the pulsar, while the extended interstellar 
medium causes the angular broadening. This model explains the otherwise- 
discordant angular and temporal effects. 

Apparent angular sizes were also measured, at 111 MHz only, for four 
other pulsars: PSR 0834+06, C950+08, 1133+16, and 1919+21. Each pulsar 

is only slightly resolved, but within the experimental uncertainty the 
apparent sizes of these pulsars and of the Crab nebula pulsar are consis- 
tent with the scattering size scaling as the square root of distance (c r 
dispersion measure). 
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CHAPTER I 


INTRODUCTION AND PURPOSE OF THESIS 

Pulsars have provided theoretical astronomers with objects which exist 
in a phy^xcal state previously mly postulated. They offer a passive astro- 
physical laboratory for the study of highly condensed matter, extreme mag- 
netic fields, and high energy plasmas. Observational radio astronomers 
have also been provided with a new phenomenon: pulsars are essentially 

point sources. But pulsars are not only point sources, they are also 
sources which very elegantly turn themselves off periodically to allow 
observers, without ever moving their telescopes, o make calibration runs 
and off-source measurements. 

Many of the spectral, temporal, and spatial features which have been 
observed in the radio radiation from pulsars are due to scattering in the 
media through which the radiation has passed. Narrow-band structure in 
the frequency domain has been studied by Rickett (19/0), Huguenin et al* 
(1969), Ewing et al. (1970), and Lang (l°71a). Temporal variations of 
the flux density (due to drifting diffraction patterns) have also been 
studied by several observers: Lang and Pickett (1970), Galt and Lyne 

(1972), Rickett and Lang (1973), and Slee et al. (1974). To investigate 
the angular structure imposed on the radiation Galt et al* (1970; and 
Vand^uberg £t al. (1973) used the technique of very-long-baseline inter- 
ferometry (VLBI) • The method of Interplanetary scintillations (IPS) has 
been used by Crcnyn (1970a), Zei3sig and Lovelace (1972), and Armstrong 
£t al* (1973) to study the angular structure. Also, the Crab nebula pulsar 
was studied by many observers, e.g. Bell end Hewish (1967), via IPS long 
before it was known to be a pulsar. 
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The project described in this thesis exploited the VLBI technique 
to achieve the high spatial resolution (~0.1 arcsec) necessary to study 
the angular structure in the radiation from pulsars. This resolution is 
obtainable only with VLBI and IPS and has added another dimension to the 
understanding of pulsars and the propagation media* This project represents 
the first application of the VLBI technique to pulsars after the early ob- 
servations by the Canadian VLBI group (Galt et^ al * 1970) • 

A* Interstellar scattering and interferometric observations of pulsars 
This section describes the phenomena associated with a simple thin 
screen model for interstellar scattering. The VLBI experiments described 
in the next chapters provide tests of this model's predictions. The agree- 
ment of observations with the model is discussed in Chapter V. 

With an interferometer we can measure the total time-averaged power 
from a pulsar and obtain a certain amount of correlated flux. From the 
pulsing power we may also obtain what we call a "correlated pulse profile." 
The correlated profile might be expected to appear quite different from 
that observed with a single telescope for the reasons discussed below. 

1, Scattering . The most important effect on the pulse profile is 
the fact of scattering of the pulsar's radiation in the interstellar medium. 
Assume for simplicity that the density irregularities which are responsi- 
ble for the scattering are distributed uniformly along the line of sight 
to the pulsar. For many purposes we can then represent the interstellar 
medium simply as a in phase-changing screen (figure 1-1) . Two effects 
of the scattering should be . *i< 3t, the radiation from the pulsar 

will be received fro. 1 s ' * - , , J.ar extent 6 which is much larger than 

s 

the intrinsic size of Second, since the radiation which is 
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received from the edges of the 
scattering cone has travelled a 
longer path than the radiation 
which comes straight through 
unscattered , a narrow pulse 
emitted at the pulsar (as in 



Figure 1-1. Thin screen 
scattering configuration. 


figure I-2a) will be broader when 

it is received at the earth. As the time delay, x, approaches (figure I- 2b) 
and exceeds (figure I-2c) the pulse period, P, the pulse shape broadens and 
a baseline of continuous flux 


is added to the profile. 
Eventually an average con- 
tinuum value of the flux is 
reached. This phenomenon was 
first proposed by Cronyn (1970b) 


i , t<* P 

I -P-l 


. T»P 


Figure 1-2. Schematic pulse 

as an explanation for the ap- trains for various time delays, 

parent decrease in the amount of pulsing flux which was observed from the 
Crab nebula pulsar at low frequencies. The amount of broadening is 
wavelength-dependent such that at longer wavelengths the average time 
delay and scattering angle both increase. 

2. Interferometric observations . Now if we observe the pulsar with 
an interferometer which has a fringe spacing on the order of the scatter- 
ing size 6 , we would partially resolve the apparent source and could then 
s 

estimate its angular size. Cohen and Cronyn (1974) point out that in order 
to observe an apparent angular size of 0 g rather than <|>» we must also use 
a bandwidth greater than 1/t or a time constant longer than the character- 
istic time for the drifting diffraction pattern. Either of these two 
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conditions insures that the rays arriving from various directions within 
the scattering cone are incoherent. 

The pulse profile which the interferometer observes would originate 
from the central area of the scattering disk where the radiation has not 
travelled as long a distance as in the outer parts. Therefore the pulse 
profile observed with an interferometer which partially resolves the scat- 
tering disk should be narrower than that observed with a single telescope. 

Another way of considering this phenomenon is to note that the appar- 
ent angular size of the source increases as the pulse progresses. At the 
start of a pulse, the radiation comes from a point source; at a later phase 
in the pulse period the radiation arrives from an annulus around the center; 
finally near the end of the period the radiation comes from near the edges 
of the scattering cone and we observe a large source. This changing appar- 
ent source size means that the fringe visibility, y, should decrease as a 
function of pulse phase, and that the pulse from a partially resolved 
source should appear to be narrower than that from an unresolved source. 

Figure 1-3 is a schematic drawing of an unresolved broad pulse and a 


partially resolved pulse (a) , 
and the fringe visibility as 
a function of pulse phase (b). 
Using an interferometer with 
very high resolution and high 
sensitivity the unscattered 
pulse shape could presumably 
be observed. The visibility 
function for pulsar radiation 



Figure 1-3. Pulse shapes (a) 
and fringe visibility func- 
tion (b). 


is discussed in more detail by 


Cronyn (1974). 



5 


Another cause for an expected difference between the corrrl<ned pulse 
profile and the single-telescope profile is due to the pu^se baseline or 
zero level. Single-telescope observers usually choose a baseline Level 
at the minimum of the pulse profile, thus assuming uhat the pulsar actu- 
ally "turns off" between pulses. If the pulse profile has both pulsing 
and non-pulsing power (which may be due to scattering, as described above, 
or to the intrinsic radiation pattern of the pulsar) then this method ig- 
nores the continuous source. With single-telescope observations, a pro- 
cedure using on- and off-source measurements would be required to detect 
the continuous component. However such an observation would requir ^ hi fc h 
sensitivity because most pulsars are weak sources and any underlying con- 
tinuous component would be still weaker and more difficult to detect. 

With an interferometer, however, the large-scale structure of the sky 
background Is washed out. The pulsar has a small angular size and its 
radiation, including both non-pulsing and pulsing components, will be 
correlated over a long baseline. Hence the level of zero correlated 
power is the true zero level of the pulsar f s power* 

In the case of the Crab nebula pulsar, the nebula makes single- 
telescope observations even more difficult, A long-baseline Interfer- 
ometer has the unique advantage of suppressing the strong radiation from 
the spatially large nebula and giving the ability to observe only the small 
compact source which is the pulsar. 

B. Summary of thesis and thesis work 

This thesis is concerned with low-frequency (100- to 200-MHz) VLBI 
observations of several strong pulsars, concentrating on the Crab nebula 
pulsar. The observed angular structure is interpreted in terms of a model 
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in which toe scattering of the Crab nebula pulsar's radiation occurs in 
both the nebula and the interstellar medium. The experimental techniques , 
from data-taking through processing to calibration, are described i.a chap- 
ter II; attention is focused on the peculiarities presented by * -tsars as 
sources observed with an interferometer. The next two chapters discuss 
the results of the Crab nebula pulsar observations: first the average 

total flux (Chapter 111) , and then the pulse profiles and pulsing power 
(Chapter IV) . Each of the various models of the Crab nebula and the in- 
terstellar medium which have been proposed to explain the observed phenom- 
ena is discussed in Chapter V. Comparisons are made between predictions 
of the models and the new results presented in this thesis. Chapter VI 
presents the results of observations of several oth-ir pulsars. Finally 
Chapter VII summarizes the results and presents suggestions for further 
work. 


The series of experiments on which this thesis is based involved a 
multitude of people: figuratively, if not literally, a cast of thousands. 

This joint project involved people at the University of Maryland, Goddard 
Space Flight Center, Arecibo Observatory, and the Naval Research Laboratory. 
It is difficult to separate the tasks of individuals in a VLBI group, but 
I shall attempt to outline the areas in which I played an active role. 

Then, in the interest of coherence, the remainder of the thesis will not 
distinguish my separate work. 

The original experiment was conceived by T. A. Clark, W. C. Erickson, 
and J. J. Broderick, and I entered the project after the first observations 
had been analyzed (Erickon et al. 1972). From then on I was involved in 
every aspect of the planning, observing, processing, and analysis operations. 



7 


During the thesis work, my "official" advisor was W. C. Erickson; my work was 
also directed by T. A. Clark who acted as my "unofficial" advisor. 

During the course of the experiments a wide variety of objects was 
observed: pulsars, compact galaxies and quasars, and supernova remnants. 

The responsibility for detailed data analysis fell to one person for each 
of the three source categories. The observations of extragalactic objects 
were studied by G. M. Resch and the results of the analysis are presented 
in his thesis (Resch 1974). Observations of the supernova remnant Cassi- 
opeia A were analyzed by L. K. Hutton (Hutton et al. 1974). The Crab nebula 
pulsar and the other pulsar observations were analyzed by me and the results 
are contained in this thesis. 

For pulsar data 1 made extensive modifications to the main processing 
program (VLBI1) to include the ability to gate the correlation routine. 

This also required that the normalization scheme of the second processing 
program (VLBI2) be completely revamped which I also accomplished. I also 
made many additional program innovations which were not directly involved 
with pulsar observations because I had assumed responsibility for maintain- 
ing the software system for all of the VLB1 data processing done at Goddard* 

As an aid to scheduling observations for VLBI experiments, I wrote the 
PRESKED and SEED programs which are briefly described in Appendix A. These 
two programs are now used by several other VLBI groups. After observing 
priorities were agreed upon in group discussion, I did the major part of 
all the detailed scheduling for the experiments discussed in this thesis. 

I participated in setting up the equipment for every experiment, and, during 
the course of the seven observing sessions, I assisted in tape-hanging at 
each of the three telescopes. 
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Once the data tapes were returned to Goddard after the recording ses- 
sions, I took charge of the data processing operation for all but one of 
the experiments. This work involved sorting and pairing all ot the data 
tapes, preparing jobs for the computers, and finally examining the output 
of VLB II and VLBI2 to decide if re-processing was necessary. For each 
experiment I also processed several runs to search for the proper clock 
settings. Normalization of the raw correlation coefficients and flux 
calibration schemes were agreed upon by the group in meetings and discus- 
sions in which I participated. Application of these schemes to the data 
was done by G. M. Resell and forms a part of his thesis. 

Analysis of the pulsar data was done entirely by me with suggestions, 
advice, and discussion by the group. This thesis represents my study of 
the pulsar data. 

C. Other research 

Before m> thesis work was begun, I participated in other VLBI projects 
including the many "Oktoberf est 11 relativity experiments and several other 
experiments in what is now called the "Quasar Patrol" project (Knight et al. 
1971; Whitney et al.1971). 

During the time that the low-frequency experiments were being done, 

I began a series of VLBI experiments at 2.3 GHz (13 cm) in order to measure 
accurate positions of pulsars relative to extragalactic sources. The col- 
laborators were G. S. Downs and P. E. Raichley of the Jet Propulsion Lab- 
oratory, and T. A. Clark, W. C. Ericsson, and myself. One purpose of the 
program was to compare positions measured with comparable accuracy using 
VLBI and pulse timing techniques in an attempt to resolve some discrepant 
measurements. Another purpose was to examine the positions derived from 
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experiments spaced several months apart for evidence of proper motion and/or 
parallax of pulsars. Manchester et_ al^ (1974) have observed proper motion 
of PSR 1133+16 using pulse timing techniques. 

Two observing sessions were held in 1972 using the Goldstone antenna 
and the NRAO 42-m telescope. The data were reduced and analyzed to deter- 
mine fringe visibilities of various sources. The Vela pulsar, PSR 0833-45, 
is found to be partially resolved. Its measured angular size is presumably 
due to scattering along the line of sight through the Gum nebula and the 
rest of the interstellar medium. More accurate normalization factors are 
being derived to improve the estimate of fringe visibility. The data will 
also be analyzed to yield relative position measurements, which can hope- 
fully be combined with the results of other interferometric studies of 
pulsars now in progress at CalTecf md NRAO. 



CHAPTER II 


EXPERIMENTAL METHODS 
A. Description of experiments 

Since 1970 there have been several meter-wavelength VLBI experiments 
using a variety of baselines in the United States. Clark and Erickson 
(1973) have summarized some of these experiments. Figure II-l is a map 
showing the baselines for the 
experiments discussed in this 
thesis and table II-l lists 
the dates, stations, obser- 
ving frequencies, and reso- 
lutions . 

1. Earliest experiments . 

One early meter-wavelength VLBI 
experiment (in 1970 January at 
121 MHz from Maryland Point to NRAO) provided the first low-frequency VLBI 
data on the Crab nebula pulsar. The technique for obtaining the Fourier 
components of a pulse shape was developed by W. C. Erickson and T. A. Clark 
and was first used for these observations. The data have been analyzed and 
published by Erickson et al. (1972). 

Another meter-wavelength experiment (1971 March at 144 MHz between 
Owens Valley, Vermilion River, and Sugar Grove) used one partly-steerable 
and two fully-steera’ *e telescopes which enabled extensive hour-angle cov- 
erage of several sources. Figure II-2 shows the (u,v) -plane coverage (for 
the Crab nebula) of these three baselines along with those of later experi- 
ments. The pulsar's scattering disk was substantially resolved with the 



Figure II-l. U.S. map showing 
baselines. Refer to table II-l 
to identify stations. 
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Low-frequency VLBI experiments 
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longest baseline and the range of fringe visibilities gave some indication 
of the brightness distribution. The experiment demonstrated that inter- 
stellar scattering does occur and that meter-wavelength VLBI is technically 
feasible. However, since all sources except the Crab nebula pulsar were 
partially resolved with the longest baseline it was also shown that more 
sensitivity was necessary in future projects* 

2. Latest experiments . With the results of the early experiments 
in mind, a multi-purpose meter-wavelength VLBI program was begun in late 

1971 and spanned one and one-half years (see table II-l). There were five 
two- or three-day observing sessions monthly from 1971 November through 

1972 March, and week-long sessions in 1971 December and 1973 February. 

The telescopes used were Arecibo (305 m) , NRAO (92 m), and Sugar Grove 
(46 m) . They form two long baselines (2550 km) and one short baseline 
(50 km). 

a. Primary purpose of experiments . The original purpose of 
the monthly experiments was simply to monitor the apparent angular size 
of the Crab nebula pulsar and compare any observed changes to the varia- 
tions in the pulse shape which were being observed at Arecibo. These 
pulse shape variations have been reported by Rankin and Counselman (1973). 
Comparisons between the size variations and the different pulse shapes 
could yield information on the location and dynamics of the scattering 
medium. The VLBI experiments were done simultaneously with the regular 
Crab nebula pulsar monitoring observations at Arecibo. 

Another purpose of the experiments was to determine the wave- 
length dependence of interstellar scattering by observing several strong 
pulsars, especially the Crab nebula pulsar, at two frequencies, 196 and 
111 MHz. In the later experiments, the observing frequency of 74 MHz was 
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Figure II-2. (u,v)-plane tracks for Crab nebula. Each track is labelled to 
Indicate requency and baseline. Symbols for each baseline are contained in 
table II-l. 
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added and 196 MHz was deleted because it turned out that interstellar 
scattering effects were not strong enough at the higher frequency. Various 
other aspects of the Crab nebula pulsar were also investigated. Enough 
Fourier cemponents were detected to allow reconstruction of pulse shapes. 
Also of interest was the polarization of the pulsar at low frequencies 

The Crab nebula pulsar was the primary source in all of the ex- 
periments, and in order to minimize ionospheric and interplanetary prop- 
agation effects the observations were scheduled in the winter when the 
Crab nebula is near the anti-solar direction. 

b. Other sources . A large number of extragalactic sources 
suspected of having, or known to have, compact components were also part 
of the program About 120 sources were observed at least once. Several 
of the stronger sources were also used for polarization and total power 
calibration observations. The intent was to obtain low-frequency spectral 
information on the compact components alone by resolving large structures 
such as halos. For many sources, the synchrotron self-absorbed region of 
the spectrum should occur near our observing frequencies. Another purpose 
of the extragalactic source observations was to compare the VLLI-observed 
angular sizes to those derived from interplanetary scintillation (IPS) 
techniques. The extragalactic sources will not be discussed further in 
this thesis except as calibration sources. The results of the observa- 
tions are presented by Resch (1974) i 

In the 1972 December experiment, a selection of supernova remnants 
was observed in a search for compact sources such as the pulsar in the Crab 
nebula. No such sources were found, although fine structure similar to that 
seen at centimeter wavelengths was found in Cassiopeia A. The results of 
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these observations are reported by Hutton et al . (1974), and are not dis- 
cussed in this thesis. 

c. Scheduling considerations . Due U the limited hou.:-angle 
coverage of the moveable feeds of the NBAO and Arecibo teles .opes, the 
amount of time a source was visible at all three stations was never more 
than about >_ne hour. Figure 11-2 showed the pitifully short (u,v) track* 
for the Crab nebula cn these three baselinas. The primary sources (;' > 

Crab nebula pulsar, five strong pulsars, and four calibration sources) 
were observed during their entire period of mutual visibility every day. 

For the remainder of the extragalactic sources, 't least two runs were 
scheduled on each of two different days. The redundancy was irsurance 
against sporadic interference, tapes with excessive errors, and '1 the 
possible problems which plague o.. . in the heat of battle. Sched- 

uling programs and other preparation procedures are described in Appendix A. 

There is as much data from each of the last two observing sessions 
as from the first five sessions combined. Also, changes in the character- 
istics of the Crab nebula pulsar did not occur during the first winter. 
Normalization and calibration procedures were improved between the two 
winters also. For these reasons, differences will often be noted between 
results of the two "sets" of experiment 0 where the first set comprises the 
first five sessions in winter 1971-2 and the second set comprises the two 
sessions in winter 1972-3. 

3. Recording system . For all of the low-frequency experiments the 
Mark I VLBI recording syBtem was used. It has a 350-kHz bandwidth and 
records a clipped and sampled video signal on standard computer tapes; 
one tape holds three minutes of data. The data can be processed with 
either a general purpose computer or a special purpose processor such 
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as the one at Haystack Observatory. The system Is described more completely 
in Appendix A. Another recording system is use^ by many other VLBI groups. 

It is known as the Mark II system (Clark 1973) and has a 2-MRz bandwidth. 

The data is clipped and sampled and recorded on video tape, ne tape lasts 
from one to five hours. The data must be processed on special purpose pro- 
cessors. 

There were several reasons for choosing the Mark I system over the 
Mark II system. The two major motivations were 1) the certainty that ac- 
curate correlation amplitudes could be extracted from the raw data, and 
2) the sophisticated software system available for data reduction. The 
accuracy of the programs was recently tested by recording a partially- 
correlated "man-made” signal on two recorders at the same site. Norma], 
process ng of the data produced correlation amplitudes within 5Z of the 
known fraction of correlated signal. This assured us that the recording 
and processing operations do indeed enable us to recover the true signal. 

The software system of the MIT/Goddard reduction programs has been 
written with exhausting attention to detail (Whitney 1974; Hinteregger 
1972). In particular, the phase is carefully tracked throughout the re- 
duction so that the recovery of pulsar Fourier components is possible. 

The ability to activate the correlation routine only when a pulsar was 
"ca" was added to the programs in order to obtain increased sensitivity 
on pulsar observations. None of these features was availaole with the Mark II 
system at the time this project vas begun, and in fact several approxi- 
mations which make it impossible to do these types of data reduction were 
built into the Mark If pre-processor. 

The other reason for using Mark I was that the integration time and 
bandwidth advantages of the Mark II system could not be exploited. At 
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frequencies between aboat 100 and 200 MHz tnere are numerous bands with 
interfering signals from radio and television stations. The Mark I band- 
width is just narrow enough to fit between most of the interference, but 
we probably could not have used the entire bandwidth of the Mark II system 
since this requires a clear band simultaneously at all the stations. Also 
at low frequencies phase fluctuations due to the ionosphere and interplan- 
etary medium limit the coherent integration time to a few hundred seconds 
or less so that we could not have made longer coherent observations with the 
Mark II system. 

B, Description of processing 

1. Program VLBI1 . The data tapes from all the experiments described 
in the previous section were carried or mailed to Goddard Space Flight Cen- 
ter and processed on the IBM 360/91 and 360/75 computers using the programs 
called VLBI1 and VLBI2. The function of each component of the interferom- 
eter and of the programs is shown in figure II-3. Details concerning the 
programs, the algorithms, approximations, and the theory of processing can 
be found in the theses of Whitney (1974) and Hinteregger (1972). The ideas 
and equations necessary to understand the basic processing and the Fourier 
component method are considered next. 

a. Ephemeris dot a . A set of ephemeris data peculiar to the 
experiment being processed is provided to program VLBI1 (figure II-3 at 
location a) . The data set contains such information as the position of 
each source observed, the locations of the telescopes, the total local 
oscillator frequencies, any instrumental time delays, and the individual 
stations’ clock epochs. At low frequencies the effect of ionospheric 
refraction is often important, and it acts as an additional, but unknown, 
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Figure II-3. Block diagram of functions of recording systems 
and processing programs. Locations labelled (a) through (f) 
are referred to in the text of sections B-l and B-2. 
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time delay. The a priori values of the parameters in the data set are as 
accurate as possible, but small errors in these values or in instrumental 
or propagation effects are not critical because the final program searches 
for the maximum correlation over a range of values. 

b. Phase . It is most important to understand vhat is meant by 

phase throughout the processing programs: how is the a priori phase cal- 

culated, how is the observed phase determined, and what information does 
the residual phase contain? 

The expression which is used to calculate the ja priori phase is 
derived in Appendix B and we will use the result: 

(<t U ‘ *. * f4, « + % “lc • U 1 - 1 * 

The first term, <j> Q , is called the "fundamental phase" which is due to the 
geometrical delay. The second term accounts for the phase due to the rate 
of change of the geometrical delay as the earth rotates; the total fringe 
rate is f and the phase is counted during the time elapsed, At^, since a 
reference time t Q . The last term includes any frequency offset between 
the two local oscillators, Au^; the phase is counted during the time 
elapsed, At^, since the start of the oscillator at t^. 

c. Correlation process . The observed correlation function, 
which describes the fringes, is obtained by multiplying the bit streams 
from two data tapes (figure II-3, at c) after shifting by the number of 
bits corresponding to the total time delay (figure II-3, at b). Because 
the time delay "resolution function" has a width ~(2ir * bandwidth) \ a 
range of ±3 bit delays centered on the nominal delay is used to define 
the correlation function. Then the seven-point function is mixed with 
"fake" fringes (figure II-3, at d) having the _a priori calculated phase 
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of equation II-l. This process is called M fast fringe rotation" and is 
done to reduce the fringe frequency and time delay to small values which 
are more easily handled. 

The phase, after mixing, is: 


4> - 

obs calc 


f R 4t R + “ T R 


(II-2) 


where f R = residual fringe rate, 

t r = residual time delay, 

u) = total local oscillator frequency. 

The quantity wt is also called the residual phase, $ , and it is made 
K K 

up of two terms: 


= A<f> + <J>. „ (II-3) 

R T o mstr 

where A4>^ = resicual fundamental phase, 

4>. = instrumental phase, 

instr r 

The instrumental phase contains terms such as clock errors and instrumental 
time delays which change slowly or not at all with time. The unknown iono- 
spheric delays also enter in this term. The instrumental phase is indeter- 
minate since it also includes the absolute phase difference between the two 
local oscillators; it can be eliminated only by using differential measure- 
ments • 

The correlation function is next averaged (figure II-3, at e) 
over 0.2 sec (one record) which restricts the residual fringe rate to 
±2.5 Hz. The resulting complex correlation functions for the nominally 
900 records in one run are written on an output tape (figure II-3, at f) 
along with all the ephemeris information relevant to the observation and 
many numbers calculated by the program. The fringe amplitudes and phases 
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will be manipulated and averaged by the next processing program to obtain 
the best estimates of the fringe rate, time delay, fringe phase, and fringe 
amplitude. 

2. Program VLBI2 . a. Residual parameters . The second processing 
program is called VLBI2. It reads the -~900 complex correlation functions 
into an array which represents the complex observable space. The data 
array is Fourier transformed irto the (residual time-delay, residual fringe- 
rate) domain; this array extends ±4.5 usee in the time-delay dimension and 
±2*5 Hz in the frequency dimension. The peak of the correlation function 
is searched for in this plane, and the maximum is found at some value of 
x and f . The best-fit values of t and f are determined by interpolation 
using a finely-spaced Fourier transform around the preliminary value* 

b* Fringe amplitude and phase * The seven points in each cor- 
relation function are then averaged using the best-fitting value for t 

K 

to adjust the phase across the bandpass* The resulting complex numbers 
are the amplitudes and phases of the fringes for each record. The shortest 
coherent integration time is one record (0.2 sec) and averages of fringe 
amplitude and phase for various longer coherent integration times are also 
calculated by the program. The choice of an appropriate integration time 
is discussed in part C-3 of this chapter (page 40). 

c* Uncertainties * The maximum fringe amplitude is not a single 
point but rather an "island" surrounded by smaller islands which are the 
noise and sidelobes of the response pattern* Figure II-4 is a schematic 
drawing of this island. The shape of the pattern is approximately a sine 
function and the width, i*e* the extent of the island, is an indication 
of how well the fringe amplitude and phase can be determined. For one 
observation, the fringe rate can be estimated to an accuracy of at least 
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(2it x T sec)”* where T is the 
coherent integration time. For 
a 180-sec average, the uncer- 
tainty in fringe rate is 
~1 m Hz . The t ime delay 
can be determined to at 
least (2ir x 350 kHz) -1 
~0.5 ysec. These esti- 
mates can be improved if the signal to noise ratio is very high. 

3. Pulsars: Fourier component method . This section discusses how 

the response of the interferometer changes when the source is a pulsar 
which has a characteristic time signature. Then the method used to re- 
cover the Fourier components of the pulse shape is described. The approach 
used in this section is different from that presented in the paper by 
Erickson et al . (1972) . 

a. Pulsar intensity. We can specify the pulsar's intensity 
as a function of time by 


f R 



Figure 11-4. Correlation peak 
in the (* R »f R ) plane. 


^(t) - I g + IjCt) (II-4) 

where l^,(t) * total intensity, 

Ig * steady or non-pulsing component of the total intensity, 
Ip(t) ■ intensity in the pulse train. 

We can further represent the total Intensity as a Fourier sum, 

00 

* T (t) ■ A n exp 2iri(nFAt R + $ n ) 


(H-5) 
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where = amplitude of the nth Fourier component of the pulse shape 

[note that A q Includes Ig and is the time average of 1^(0], 
* phase of the nth Fourier component, 

F * 1/P ■ pulse repetition rate, 

Atg ■ elapsed time since the reference time, t Q . 

The various components of the total intensity are illustrated in figure 


II-5. 

This is not a complete 
description of the source 
radiation however, because 
all pulsar radiation is dis- 
persed in its passage through 



Figure II-5. Components of the 
total intensity of a pulsar. 


the interstellar medium. The 

effect of dispersion is to shift the pulse later in time by an amount de- 
pending on the frequency and in such a way that pulses arrive later at 
lower frequencies. To illustrate this, figure II-6 shows schematically 


some pulses in the (time, fre- 
quency) plane. The pulses 
sweep across a given pass- 
band at a rate Au/At, which, 
over a small range of fre- 
quency, can be approximated 
as a constant. 



Figure 11-6. Dispersed pulses 
in the (time, frequency) plane. 


Au/At 


-u J -D = -F/G 


(II-6) 



where D - dispersion constant, measured as the difference in the time of 
arrival of the same pulse at two frequencies, 

3 

1/G = spacing of the pulses in the frequency domain (G = 2FD/w ), 

1/F « spacing of the pulses in the time domain (pulse period) . 

Since the entire pulse is delayed, the pulse shape is not changed and each 
Fourier component is affected the same way. The time delay due to dis- 
persion is 

t D = DA (l/o) 2 ) (II-7) 

2 

which can be represented as a phase shift of ~2TTinFD/u> . The expression 
for the total flux including the effects of dispersion is 

I T (t) = £ A n exp 2iri(nFAt R + nGw + <j> n ) . (II-8) 

n=- 

b. Response of interferometer . The response of an interferom- 
eter is the product of the source intensity and the fringe pattern, 

R(t) = J^Ig + Ip(t) j k exp 2iri4>(t) . 

When <j>(t) is given by equations II-2 and II-3 for the residual phase, we 
have 

R(t) = £ A n exp 2iri [(A(fi o - nGw) + (f R - nF)At R + ^ + 4> n 1 • 

n— oo L J 

(H-9) 

Notice that the residual fringe rate is given by (f D - nF) . This means 
that for a pulsar there are many sets of fringes, one for each value of 
n, in contrast to only one peak for a steady source. The fringes are 
separated from each other in fringe itequency by the pulse repetition rate 
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and their amplitudes are modulated by the amplitude of the nth Fourier 
component, A r . Also, since the pulsar is dispersed, there is a phase 
change across the passband. This is shown in equation II-9 by the fact 
that the residual fundamental phase is given by (A<J> q - nGw) . In the 
transform domain a phase shift corresponds to a time delay, so that each 
stc of fringes occurs at a different delay. Figure II-7 shows the (? D ,f ) 
plane for a dispersed pulsar 
having F « 1 Hz and G * 1 ysec. 

Appendix C shows pictorially 
the various sets of fringes. 

c. Recovery of A^ and <f> . 

To obtain the amplitude and phase of 
the nth Fourier component, a time 
delay of nG and a fringe rate off- 
set of nF are inserted during pro- 
cessing so that the main response of the interferometer is centered on the 
nth set of fringes in the (t^, f R ) plane. Since one processing pass through 
VLBI1 covers an area in the (t , f ) plane of about 9 Usee by 5 Hz, all the 
fringe islands which occur within this domain can be found by searching with 
program VLBI2 at the appropriate locations. Additional processing passes 
are required for fast and/or highly dispersed pulsars and for higher order 
components (see the pictures in Appendix C) . 

In practice, obtaining A^ and is done in the following way. 
Returning to equation II-l for the calculated phase, we set ■ nF, 

insert an instrumental time delay * nG, and specify t^ Q as some UT near 
the start of the tape. (The value of t^ 0 actually depends on the exact 



Figure II-7. Multiple sets 
of fringes for a dispersed 
pulsar. 



manner in which time is referenced in the programs.) The calculated phase 
we obtain is then: 
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W *0 + £it R + ” F4t L0 • 

The observed phase of the nth Fourier component is: 

*(t) - nF(t - t ) (11-11) 

u n 

where t is the time at which the sine wave of frequency nF has its maxi- 
mum or epoch. The output of VLBI2 is the residual phase of the nth compo- 
nent: 

- A<J> + * + nF(t - t T .) , (II-12) 

nR o T instr n LO' * 

and for n * 0 the residual phase is given by equation 11-3. 

Differential measurements between the residual phases for n ■ 0 
and n f 0 allow the instrumental phase and the residual fundamental phase 
to be eliminated. The "differential residual phase" of the nth component 
is 




nR 


OR 


‘ nF(t n ' ^ 


(11-13) 


d. Pulse shape reconstruction . To reconstruct the pulse shape 
the differential phases of equation 11-13 and the corresponding amplitudes 
are added in a Fourier sum. Since t^Q is known in UT, the residual phases 
contain information on the epoch of each component with respect to UT, and 
the reconstructed pulse contains information which can yield the pulse 


time-of-arrival in UT 
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The Fourier components are summed according to the equation 

00 

I T (t) » £ A n exp 2TTi(A<J> nR - nFt) . (11-14) 

n»- oo 

In practice, the number of Fourier components which were detected in the 
processing Is finite so the sum is actually from -N to +N where N is vari- 
able but always less than 10, The function I^(t) is interpolated so that 
the resulting pulse shape is smooth, although it must be remembered that 
there are only N + 1 independent points in the pulse profile. For the 
computer program which reconstructs and plots a 100-point pulse, the sum 
is represented as 

+N 

VV ■ £ A n exp 2iri(A<k - nt^/100) (11-15) 

1 n*-N 

where i = 0 to 99. 

The pulse shape which results from this process is an average 
pulse shape unsmoothed by dispersion effects. This profile is equilvant 
to that which would result from single- telescope observations of the aver- 
age pulse shape over three minutes having no dispersion and having the 
sensitivity of a 350-kHz bandwidth system. 

e. Zero-baseline interferometer . Because there is frequency 
structure within the Mark 1 passband, we can also obtain information about 
the pulsing power by auto-correlation. The data is processed in the same 
way as for cross-correlation in that time-delay and fringe-rate offsets 
are cr.-.'lied, but only one tape is used. The Fourier components n * 1, 2, 

... H are obtained. Because this is an auto-correlation process the n ■ 0 
component has a correlation coefficient of unity and the amplitudes of the 
+n arid -n components are identical and the phases have opposite signs. The 
pulse shape is reconstructed in the same manner as for the cross- correlation 



profiles, the only exception being that an arbitrary-level n = 0 component 
must be added to raise the baseline of the pulse* The components obtained 
this way and the resulting profile are called "zero-base Line 11 romponents 
and profiles* Usually the Arecibo Mark I data tape was used for the zero- 
baseline processing because the sensitivity was the highest at this station* 

4* Pulsars: gating method * a* Description of processing * From 

pulsar VLBI data we can obtain either the Fourier components, as just de- 
scribed, or the time-averaged power* The interval o^er which the average 
power is measured can be restricted by "gating" the correlation routine so 
that it operates only when the pulse is "on." This makes the process more 
efficient since most pulsars are "off" more than 90% of the time. 

The gate used with the VLBI1 correlation routine is not a true box- 
car function but is actually "fuzzy." This is because 2 msec is the short- 
est segment of data for which the correlation routine is called, and if a 
pulse begins anywhere within that segment the correlation is done; simi- 
larly, the trailing edge of the gate is only defined to within 2 msec. 

The record length, 0*2 sec, is not conmensurate with any pulse period so 
that some pulses may straddle two records with the resulting loss of about 
5 msec of data, the length of the inter-record gap. If the pulse width is 
less than 5 msec, an entire oulse may be lost in the gap. These contin- 
gencies are handled by the subroutine which opens and closes the gate 
The observed correlation function is normalized by the actual amount of 
data correlated so the duty cycle of the pulsar is automatically accourted 
for in the correlation coefficient. The remainder of the processing is 
the same as for a continuum source. 

b* Epoch measurement * The time-of-arrival (epoch) of one pulse 


must be known to an accuracy of about one millisecond so that the gate can 
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be placed properly over the pulse. Although it is possible to search for 
the correct epoch by placing the gate at different trial epochs until the 
maximum correlation is found, this is time-consuming and impractical. 

To monitor the average pulse shape during the observing sessions, 
the IF signal was sent to a multi-channel signal averager which wes trig- 
gered by a clock running at the pulse repetition frequency. The stored 
signal was displayed on an oscilloscope and a polaroid picture was taken 
to record the profile. The time at which the clock had bean staged was 
recorded so that the hi of the start of the trace was known. The epoch 
of the pulse could be read off the picture in units of percentage of the 
pulse period, giving the pulse epoch relat .ve to UT. This time was re- 
corded and later used in the processing program to calculate the time at 
which the gate should be opened. 

Besides knowing the epoch of the pulse, another requirement for 
using the gating method is thit there be only one pulse in the Mark 1 
passband at any one time. This is the case for all of the It w-dispersicn 
pulsars we observed. However, the gating method as described above cannot 
be used for the Crab nebula pulsar because its dispersion measure is so 
high and period so short that several pulses (e.g. at 111 MHz) are 
present across the band simultaneously, il we had i •ned the g a ’ng method 
for this pulsar, we woulc add up "on" and "off" signals when we integrate 
across the passbard, which is what the method tries to avoid in the first 
place. For the Crab nebula pulsar, the pulsing power is recovered by de- 
dispersing the signal as described in the previous section. 

c. Sensitivity and flux measuremen t** > The intensity of pulsar 
radiation is expressed in a variety of ways. Pulsar spectre are usually 
presented as plots c log E vs. log w where E is the " average pulse energy 11 
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-2 -1 

in units ot J m Hz This quantity is not dependent on the details of 
the pulse shape, it is simply a measure of the average received energy. 

The intensity is also given in terms of the “peak flux density,” = E/^, 
where w^ is the equivalent pulse profile width. Yet another quantity is 
the “equivalent continuum flux density,” * E/p, where P is the pulse 
period. is the total time-averaged flux density, that is, the flux 
which the pulsar would have if it were smeared into a non-pulsing contin- 
uum source. By using the gating method we measure the flux density aver- 
aged over the gate width: * E/w^; nominally we set > w^ to include 

most of the pilse profile. 

Gating increases the signal-to-noise ratio in the following way. 

If we do not gate, then we must detect a continuum source signal in the 
presence of a noise level With a gate, the integration time is de- 

creased by a factor w /P, so the noise level is increased to a « a /P/v , 

g peg* 

tut now the signal to be detected is larger, - S c ?/w^. Therefore, the 

signal-to-noise ratio is enhanced by a factor /P/w if gating is used. 

g 

Typically, P/w £ 20, so the signal-to-noise ratio improvement i > at least 
g 

a factor of 4.5 over the case where pulsars are treated as continuum sources. 

We can recover the correlation c<; efficient of a single pulse pro- 

vided that the pulse period is less than 0.2 sec, the shortest coherent in- 
tegration time. The detection limit for single pulses is very high however 
because the integration time xs only w • 

s 

d. Replication of frlnRes . In the previous section on the 

Fourier component method It was noted that for a pulsar there are many 

sets of fringes spaced by the pulse repetition frequency. For the Crab 
nebula pulsar F ~ 30 Hz so the fringes are well-separated from the ambi- 


guities which occur due to sampling by the Mark I recording system. However , 
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for pulsars with F — 1 Hz the different sets of fringes are closely spaced 
in the fringe rate spectrum and they could be mistaken for the zero-order 
fringes if not properly taken into account. The search by program VLBI2 for 
fringes is restricted to the range ±F/2 so that these ambiguities are 
avoided. Appendix C contains pictorial representations of a sample of 
data, and the corresponding fringe rate spectrum, to illustrate the rep- 
lication. 


C. Calibration procedures and correction factors 

1. Flux scale . a. General method . All of the data from the 
seven experiments were normalized by the following general method. Since 
the correlation coefficient, o, is a fraction of the mean total system 
temperature it car be expiessei as a temperature by multiplying by the 
mean temperature: 


T - p/f 7 
corr v sysl 


sys2 


(11-16) 


Measurements of each antenna's efficiency, e, expressed in Kelvins per 
Jansky, can then be used to convert the temperature to flux density: 


S - T A/ e.e_ (II-17) 

corr corr v 12 

For the 1971 March experiment at 144 MHz the correlation co- 
efficients have not been converted to fluxes because little information 
on system parameters and total power records is available. However, 
because each telescope in this experiment had approximately the same 
sensitivity the correlation coefficients, in correlation units, can be 
inter compared between the three baselines. 
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The system configuration at each station varied greatly from one 
experiment to the next as improvements in receivers and observing techniques 
were gradually made. A discussion of some of the situations encountered in 
this learning process is contained in Appendix A. 

b. First set of experiments . For the first set of experiments, 
system temperatures were estimated by adding together the known contribu- 
tions to temperature: 

T 36 T (source) + T . + T . (11-18) 

sys a sky receiver 

At meter wavelengths, all but a handful of sources are weaker than the sky 
background so the source antenna temperature can usually be neglected. The 
sky background temperature was estimated by reading the contour level near 
the source position from the 178-MHz map of Turtle and Baldwin (1962) and 
scaling to our frequencies with a brightness temperature spectral index of 
-2.5. The random error in estimating the sky temperature is about 20Z. 

The error also depends on frequency because the sky background becomes 
smoother at lower frequencies and it also forms a larger fraction of the 
total system temperature. 

The receiver temperature was accurately measured at NRAO and 
Sugar Grove only in Marcn 1972. For each station, the receivers were 
assumed to be similar enough each time that the 1972 March measurements 
were valid for every other month also. At Arecibo, the receiver system 
also was changed every experiment but here the temperature was measured 
each time, and these values were used. The receiver temperature has an 
uncertainty of about 20% due to these problems. 

The total error in the system temperatures will be due to random 
errors in scaling the sky map for background temperature, and to a systematic 
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error from uncertainties in receiver temperature. Overall, the estimate 
of the total system temperature for each station is probably accurate to 
30Z. 

The dish efficiency of the 92-m NRAO telescope was carefully 
measured at 111 and 196 MHz in 1972 March (Hutton 1972), and these curves 
were adopted for all the previous experiments. The efficiency is a strong 
function of hour angle, and the normalization from correlation units to 
Jansky s will thus be dependent on the hour angle of the source at NRAO. 

At 111 and 196 *®z the efficiency is a monotonically decreasing function 
of hour angle, but at 74 MHz the efficiency is a more complicated function. 
The adopted curves for all frequencies are shewn schematically in figure 
II-8. At Arecibo, the vari- 
ation of efficiency with hour 
angle was not as severe as at 
NRAO but the effects were 
included. The efficiency 
of the telescope at Sugar 
Grove did not vary with 
distance from zenith. The 
feed systems at all three 
stations were the same each 
month so it was assumed that 
the antenna efficiencies did not vary. 

The total error in normalization from correla*-’ - i units to 
Janskys for this set of experiments is estimated to be .. :uv 207* Table 
II-2 contains typical values of receiver temperature and efficiency; 

-able A-l lists the individual values for each telescope and observing 



Figure 11-8. Hour-angle dependence 
of efficiency for NRAO 92-m telescope. 


session 
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A special problem exists for the 1972 March experiment: the 

sensitivity of the Sugar Grove station was reduced to below its usual 
value by about a factor of three. The decrease in sensitivity appears 
to vary randomly during the experiment, however, making a correction 
factor difficult to compute. For this reason only the Arecibo-NRAO base- 
line flux measurements will be presented for this one experiment, 

c. Second set of experiments . During the last set of experiments 
the system temperatures were measured for every run. This eliminates 
uncertainties due to estimating background temperatures but introduces 
reliance on a calibration noise source. The random error in one estimate 
of system temperature is about 10Z. The curves for efficiency as a func- 
tion of hour angle at NRAO and zenith angle at Arecibo were again measured 
by tracking a source across the sky and frequently moving the telescope 
beam on and off the source. This was done using Vir A at Arecibo, and 
Cas A, Tau A and Cyg A at NRAO, The three sources at NRAO gave consistent 
results • 

The error in normalization actually depends on the total system 
temperature because of the measurement method we used; a typical error is 
estimated to be about 8Z for a single observation on a weak source. The 
error is larger, about 12Z, for observations of the Crab nebula pulsar 
because the nebula greatly increases the total system temperature. Sys- 
tematic errors occur due to uncertainty in the efficiencies; these are 
estimated to be less than 10Z. Table 11-2 contains average values of 
receiver temperature for each experiment. Table 11-3 contains sensitiv- 
ities calculated using the values of table 11-2 and the equation 

Wi • 

min " 


\/ e l e 2 \/®^ 


(H-19) 
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Table 11-2 

Average System Parameters, by station 


Date 

Parameter 

196 MHz 
A N S 

Ill MHz 
A N S 

74 MHz 
A N 

S 

Winter 1971-2 

receiver temp. (K) 

350 

350 

375 

600 

170 

675 

- 

- 

- 


efficiency (K/Jy) 

1.3 

.93 

.31 

3.5 

.96 

.30 

- 

- 

- 

Winter 1972-3 

receiver temp. (K) 

- 

- 

- 






750 


efficiency (K/Jy) 

- 

- 

- 

2.6 

1.0 

\ 

.35 

4.6 

.60 

.35 


Table II-3 

System sensitivity £or weak sources, by baseline 



196 MHz 
AH AS NS 

111 MHz 
AN AS NS 

74 MHz 

AN AS NS 

Winter 1971-2 
sensitivity (Jy) 

0.5 0.8 0.9 

0.5 1.2 > .8 

- 

Winter 1972-3 
sensitivity (Jy) 

- - - 

0.5 1.3 1.8 

1.4 1.9 5.1 

Average sky 
temp . used (K) 

150 

600 

1800 












36 


where B = 350 kHz, T = 180 sec (one tape), and and T ^ were calculated 
from the receiver temperatures at each station and average sky temperatures 
given in table II-3. 

d. Flux scale . The flux scale was established by assuming 
flux density values for Cas A, Cyg A, and Tau A at our observing frequen- 
cies. Table II-4 lists the fluxes. The values were obtained by computing 
a best fit, in the least squares sense, to the available absolute measure- 
ments at many frequencies. The details of this procedure are discussed 
by Resch (1974). 


Table II-4 

Standard source fluxes 


Source 

19o MHz 

111 MHz 

74 MHz 

Cassiopeia A 

10960 Jy 

16680 Jy 

22570 Jy 

(epoch 1972.9) 




Cygnus A 

8955 

13275 

17210 

Taurus A 

1525 

1775 

1980 

Virgo A 

1030 

1575 

2145 


Systematic errors in flux measurements can only be evaluated 
through comparison with other independent observations. Each of the 
seven experiments was normalized independently and then the values of 
correlated flux were inter-compared and also compared to source surveys 
made at frequencies near our observing frequencies. The procedures we 
used are as follows. 

Since the two long baselines are of nearly equal length, the 
correlated flux of a source observed on both baselines should be the same. 
This was tested by computing the ratio of the two fluxes (AN/AS) for runs 
in 1973 February. The average ratio was less than 5% different from unity 



37 


at both frequencies. This indicates that the NRAO and Sugar Grove effi- 
ciencies have the correct ratio. Next, the observed fluxes of about six 
sources which were expected to be unresolved on the short baseline were 
compared to independent published flux measurements. It was found that 
the fluxes agreed within 10Z at both frequencies, indicating that the 
short baseline stations had been calibrated correctly. Finally, the ratios 
of the flux observed in 1972 December to that in 1973 February was computed 
and the average ratio was found to be less than 5Z different from unity. 
This shows that both experiments were calibrated the same. 

Unfortunately, we have no way of determining the systematic 
errors involved with both long baselines because no independent measure- 
ments exist. However, the agreement of our measurement of the pulsing 
flux for the Crab nebula pulsar (using zero-baseline pulse profiles) with 
the single-telescope measurements at Arecibo supports the accuracy of the 
Arecibo baselines 1 calibration. 

These arguments establish that all of the flux values presented 
here are Indeed consistent and on the same flux scale. Unknown systematic 
errors may remain and their magnitude is estimated to be <10Z. More de- 
tails on the analysis of systematic errors are given by Resch (1974). 

2. Polarization . Polarization observations with an interferometer 
are more difficult than with a single telescope because of the heightened 
Instrumental effects. Incomplete isolation of the oppositely polarized 
feeds results in an unwanted response. A right circularly polarized (RCP) 
feed responds to a fraction, f, of the left circularly polarized (LCP) 
component of the signal, and also -to a fraction, 1 - f, of the RCP com- 


ponent 
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The quantity f is called the single-telescope instrumental polariza- 
tion. For total power measurements with a single telescope it makes no 
difference that f is not zero, unless the source is polarized. However 
the response of an interferometer is reduced because the LCP and RCP com- 
ponents of the signal are not correlated. This effect can be corrected 
for if the value of f (or some other suitably calibrated quantity) for 
each antenna is known. 

When the two stations of an interferometer' observe with oppositely 
circularly polarized feeds, the interferometer responds to linearly polar- 
ized radiation of any position angle. However, due to the fact that the 
value of f for each station is not zero, the response to an unpolarized 
source is not zero. The quantity a is called the "interferometer instru- 
mental polarization" and it is the ratio of the correlation coefficient 
due to an unpolarized source as observed with similarly- and oppositely- 
polarized feeds at the two stations. The instrumental polarization is 
approximately equal to the square root of the sum of the individual in- 
strumental polarizations, f ^ and f ^ (see Appendix D). Since the values 
of f are fractions, the instrumental polarization of the interferometer 
can be considerably higher than either of the single-antenna instrumental 
polarizations. For example, to achieve an interferometer instrumental 
polarization of 15X, each antenna can have no worse them IX instrumental 
polarization. 

Values of a were obtained from observations of strong extragalactic 
sources which were assumed to be unpolarized (3C273, 3C286, 3C287, and 
3C298). Tables II-5 and II-6 give the observed values of a and the de- 
rived values of f. The instrumental polarization was very high for all 
of the experiments, and it dominated most of the attempts to measure linear 



Interferometer instrumental polarization (a) 
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pole ,. ization of pulsars. The results of the observations are discussed 
i.a the next chapters. One consequence of the large instrumental polariza- 
tion is that the correction factor for observations with similar feeds at 
both stations is of the order of five percent, comparable to other cali- 
bration effects. 

3. Noise fluctuation level and integration times . The complex cor- 
relation coefficient (fringe amplitude and phase) can be consider'd as an 
observed vector (R) which is composed of a signal vector (S) and a randomly 
oriented noise vector (N) as shown in figure II-9, where |s| - p, |n| = a, 
and | R | = p°. The noise vector 
has a Rayleigh distribution in 
amplitude and a uniform dis- 
tribution in phase; the signal 
vector is a constant. The 
resultant vector then has a 
Rice distribution which approaches 
a Gaussian distribution as the signal-to-noise ratio increases. Because we 
measure a normalized correlation coefficient with a clipped system, the 
noise is a constant and does not depend on the signal level. The rms noise 
level is defined as 

a - ■ ft /- 2 ( 11 - 20 ) 

/St 

For a full-tape coherent average, T - 180 sec, and with a bandwidth of 
B ■ 350 kHz, o » 1,4 * 10 . For a reliable detection of fringes, we 

usually require the correlation to be above the 50 level. However, a 
small residual fringe rate and time delay can be used to verify signals 



Figure II-9. Signal, noise, 
and resultant vectors. 
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as weak as about 3o. To obtain an unbiased estimate of the fringe ampli- 
tude we subtract the noise 

f °2 2l 1/2 

p - |p - o J (II- 21) 

The rms uncertainty in the phase is 

S o/p (11-22) 

for measurements with a noise level of 0 . More details and derivations 
can be found in Whitney (1974) and Rogers (1970) . 

For coherent integration times less than a full tape in length, o 
increases as the square root of the time, and the level for detecting 
correlate ">n is higher. The need for shorter coherent integration times 
is due to non-random phase noise such as oscillator jumps, or refraction 
effects. A coherent integration which included these effects would give 
a fringe amplitude which is lower than the true amplitude. 

The phase perturbations with time scales of about one minute are 
easily seen in the plots of fringe phase versus time which are output by 
VLBI2, The deviations are probably due to ionospheric "blobs" moving 
across the beam of one of the telescopes. In such cases broken coherence 
averaging should be done: the data is coherently averaged over segments 

for which the phase is fairly linear. This procedure gives a better esti- 
mate of fringe amplitude but increases the uncertainty by reducing the 
integration time. 

If incoherent averaging is employed, it is not necessary to have 
detected the signal in one short coherent integration time in order to 
estimate the correlation coefficient. Moran (1973) states that the op- 
timum estimator for incoherent averaging, in the weak signal case, is 
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the square root of the average of the squares of the coherently averaged 
quantities. The noise level for the resulting cst^jiate decreases only as 
the fourth root of the total integration time. 

In program VLBI2 the data is coherently averaged over segments of 
lengths 0.2, 2, 20, and 180 sec and then each set of averages is inco- 
herently averaged. For observations of most sources, the 20-sec averages 
were routinely used unless examination of the phase for a particular run 
indicated that either a longer or rbc ,*ter time was * ore appropriate. 

During all of the experiments we were plagued with an additional ran- 
dom source of errors excessive tape errors which resulted from a fault 
in the tape drive at Sugar Grove. Usually bad runs ce easily detectable 
because parity errors occurred on nearly every record and the fringe ampli- 
tudes were anomalously low. Such runs, when obvious, were culled. 

4. Interplanetary scintillation effects . Over all time scales in- 
terplanetary scintillations (IPS) cause phase perturbations in the wave- 
front arriving at the earth. This effect can be corrected for by using 
a model of the solar wind, such as that proposed by Cronyn (1972). He 
derived a formula for the amount of phase deviation in the wavefront, A$, 
as a function of observing wavelength, solar elongation, and coherent 

integration time. The fringe amplitude will be reduced such that only 

2 

a fraction, exp (-A<J> /2), of the true correlated power is detected, as 
if the source were partially resolved * The amount of phase shift increases 
with longer wavelength, smaller elongation, and longer coherent integration 
times. Table II-7 lists, for elongation 90°, values of Ao at 111 and 74 MHz 
for various integration times. At 180® elongation the values of A<fr would be 
smaller by a fraction of 0.7; for sources out of the ecliptic plane, A$ will , 


also be reduced 
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Table II- 7 

IPS phase deviations for e ~ 90° 


Coherent integr. time 
Records Seconds 

111 MHz 

A<Krad) exp(-A<t>^/2) 

74 MHz 

A^(rad) exp(-A<J>^/2) 

900 

180 

1.7 

0.25 

2.5 

0.05 

100 

20 

0.68 

0.79 

1.0 

0.57 

13 

2 

0.20 

0.98 

0.30 

0.95 

1 

.2 

0.03 

99 

0.05 

0.99 


For sources in the ec 1 • ptic r lane we should use not longer than 2-sec 
averages according to Cronyn’s model. Alternatively, we could correct tre 
observed correlation coefficients by the calculated factor. Short inte- 
gration times of 20 sec were routinely used, but no correction factors for 
the solar wind were applied. Large phase deviations due to interplanetary 
scintillation of 3C48 were observed when the source was at a solar elon- 
gation of 60° in 1973 February. However, the data showed that the decor- 
relation effect was not nearly as severe as predicted by the model. These 
observations are discussed by Resch (19/4). For .he Crab nebula observa- 
tions, slight decorrelation was observed with longer integration times; 
this is discussed in the next chapter. 

D. Sumnary of Chapter II 

The techniques involved in a VLBI experiment have been discussed from 
the point of view which is most useful ,t the results which are presented 
in the next two chapters. The focus in this chapter bar b-jen on acquiring 
and processing meter-wavelength pulsar data with a VLBI bystem. Part k 
described the experiments and observations; part B was concerned rith the 
processing of interferometer data to obtain estimates of correlation 
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coefficients. Finally part C presented the various calibration procedures 
adopted for these particular experiments. The data presented in the remain- 
ing chapters has been through all of the procedures outlined above s o that 
the best estimates of correlated flux have been obtained. 



CHAPTER III 


CRAB NEBULA PULSAR: TOTAL FLUX RESULTS 

This chapter presents the results of our studies of the so-called 
"compact source" in the Crab nebula (Hewish and Okoye 1965) * The spectrum 
of the compact source is very steep and its flux dominates the nebular 
emission at very low frequencies* Fluxes measured at 74, 111, and 196 MHz 
are reported here. The measurements were made with both long- and short- 
baseline interferometers so we have partially resolved the compact source 
and have made estimates of its angular size. 

In 1968 Staelin and Reifenstein (1968) discovered a pulsar in the 
Crab nebula, and since then many arguments have been presented which point 
to the identification of the pulsar with the compact source. Most of the 
discussion in this chapter implicitly assumes that the compact source and 
pulsar are identical, and some additional, and very convincing, results 
are presented to support this view. 

A. Data presentation 

1. Definitions. There are two ways to consider the signal from the 
Crab nebula pulsar: either as a periodic wayeform which can be analyzed 

4 "to its Fourier components, or as an average power. These quantities 
were shown in figure II-5 and defined in equations II-4 and II-5 (pages 22-3). 
The time-averaged pulsing power is denoted as Ip, the non-pulsing or steady 
flux as Ig, and their sum as t:he total time-averagod flux, I T . At fre- 
quencies above about 175 MHz, all of the power from the pulsar is pulsing 
so that Ig = 0 and Ip - I T . Below about 70 MHz, essentially all of the 
power Is non-pulsing, and Ip * 0 and lg * 1^. Near 100 MHz the nulsing 
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and non-pulsing components contribute about equally to the total power and 
Ip ~ I g and I p ~ 1^2. 

Considering the pulsar as a continuum source, i.e. the compact source 
in the Crab nebula, we process the data normally and obtain values for 1^. 
These results have the highest signal-to-noise ratio because all of the 
available power is being observed* The results of observations of will 
be discussed in this chapter, and the results of the observations of Ip and Ig 
are in the next chapter* 

Information on the size and spatial structure of a source can be ob- 
tained from total flux measurements* Normally, the fringe visibility of 
a source is calculated by comparing the flux observed with an interferom- 
eter to the total flux of the source observed with a single telescope 
(zero-baseline interferometer) • For the Crab nebula pulsar, the nebula 
provides a strong continuum background against which it ^s only possible, 
with a single telescope, to observe the pulsing flux. Ip. We used a short- 
baseline interferometer to resolve the spatially large nebula so that only 
flux from the compact source is correlated; in this way we obtain a measure- 
ment of I T which includes both continuous and pulsing components . The 
fringe visibility, y, is defined as the ratio of the fluxes observed with 
the long baseline and short baseline* This assumes that the source is 
totally unresolved on the short baseline. The long-baseline flux is de- 
noted by a subscript R which means partially resolved flux. We then have 

I TR " I PR + I SR Y " I TR^ I T* 

2« Data lists and plots . In table III-l are collected the one-day 
average fluxes for each frequency, the values of frinfr visibility, and 
the over-all averages. The fluxes observed with the two long baselines 
were averaged together for this table. Figure 111-1 is a plot of the flux 
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Table III-l 

Total average fluxes of the Crab nebula pulsar 
196.5 MHz 


Date 

VJy) 

_ No^ 

I TR (Jy) 

Mo. _ 

1 

Notes 


1971 Nov. 23-24 8.1±0.1 2 4.9±0.5 7 0.61±0.10 (a) 


1971 Dec. 19-20 

7.3±0.7 

6 

6.211.5 

14 

0.8510.26 


1972 Jan. 23-24 

— 

- 

7.711.2 

9 

— 

(b) 

1972 Feb. 27-28 

10.0±0.4 

4 

6.511.0 

14 

0.6510.12 


1972 Mar. 22-25 

— 

— 

6.910.5 

3 

— 

(c) 

averages 

8.311.2 
8.310.4 p.e. 

12 

6.411.4 
6.210.2 p.e. 

47 

0.7910.27 
0.7410.06 p.e. 

(e) 


Date 

VJy) 

111. 

No. 

5 MHz 

r TR (Jy ) 

No. 

y 

Notes 

1971 Nov. 23-24 

50.0H.7 

3 

31.215.8 

* 

0.6310.12 

(a) 

1971 Dec. 19-20 

45.610.8 

3 

29.413.0 

15 

0.6410.07 


1972 Jan. 23-24 

— 

- 

28.913.2 

9 

— 

Cb) 

1972 Feb. 27-28 

49.3+1.5 

4 

28.813.4 

14 

0.5810.07 


1972 Mar. 22-25 

— 

- 

23.215.1 

4 

— 

(c) 

1972 Dec. 8 

36.112.5 

13 

— 

- 

— 

(d) 

1972 Dec. 9 

34.613.8 

11 

— 

- 

— 

(d) 

1973 Feb. 25 

41.2+3.3 

10 

2' 3.8 

20 

C. 64+ 0.10 


1973 Feb. 26 

40.214.3 

11 

21 3.0 

17 

0.72+0.11 


1973 Feb. 27 

44.212.6 

11 

30.5^2.4 

23 

0.6910.07 


averages 
adjusted avg. 

40.514.6 

40.510.6 p.e. 

66 

28.913.8 
31.211.2 p.e. 

111 

0.7210.03 p.e. 

(e) 


Date 

i T Uy) 

74, 

_ Noj_ 

.0 MHz 

XjgOy) 

No. 

Y 

Notes 

1972 Dec. 5 

105.3119.4 

11 




(d) 

1972 Dec. 6 

110.4112.1 

11 

— 

- 

— 

(d) 

1972 Dec. 7 

101.6111.1 

11 

— 

- 

— 

(d) 

1973 Feb. 28 

96.013.6 

3 

30.7±2.8 

17 

0.3210.03 


1973 Mar. 2 

105.217.9 

5 

30.713.4 

12 

0.2910.04 


averages 

105.01’ 5 

41 

30.713.0 

29 

0.2910.05 


adjusted avg. 

105.01 2.1 p.e. 

33.511.7 p.e, 

» 

0.3210.02 p.e. 



Notes: 

(a) Strong Interference at Areclbo. 

(b) Sugar Grove not on the air. 

(c) Sugar Grove sensitivity low. 

(d) Areclbo not on the air. 

(e) See discussion on p. 49 . 
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Figure III-l. Total average flux of the Crab neDula pulsar. The error bars are 
the standard deviations of all the data points. Data points are listed In 
table III-l. The errors at 196 MHz are the same size as the points. 
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vs. time using the values in table III-l. The average fluxes at 196 and 
111 MHz were computed from the data listed in the tables of Appendices E 
and F. The error quoted following each flux value in the table is the rms 
deviation of the data points and should be the random error in a single 
estimate of the flux. This is the error plotted in the figure. Values 
for fringe visibility were computed using the listed average fluxes. The 
errors were calculated as the quadratic sum of the percentage error in each 
flux used in the calculation. The overall average flux and fringe visibil- 
ity at each frequency are weighted averages, and the associated error rep- 
resents the rms deviation of all of the points. The error of the "adjusted 
average" is the standard deviation of the mean value which takes into 
account the total number of data points. The "adjusted" flux is discussed 
next. 

At all frequencies, full-tape (180-sec) coherent integration times 
were used unless significant phase perturbations or a phase discontinuity 
was noted and then a shorter time was used. Since the appropriate noise 
level was subtracted from the fringe amplitudes (equation 11-21) , the cor- 
related fluxes should be nearly equal for any integration time in which 
the signal is detected. However, a calculation at 111 and 74 MHz using 
20-sec coherent integration times increased the fringe visibilities 
(because the long-baseline correlated flux was higher) by an average of 
8±4%; this is within our estimated flux error. At both frequencies the 
long-baseline fluxes are usuall}* detected at about 15 times the rms noise 
level for a 20-sec integration time. For 2-sec integration times the 
fringe visibility is increased by 19±8% at 74 MHz and 13±6Z at 111 MHz, 
but, the flux is only a the 5 o detection level for this short Integration 
time. It was decided to use the 180-sec integration times In order to 
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retain a high confidence level in the correlated flux and the calculated 
visibilities. Then, to remove the known systematic error, the average 
long-baseline flux, the fringe visibility, and the error were raised by the factor 
1.08, to yield the "adjusted averages." Because there are many runs the 
probable error is still quite small. 

This systematic effect is indicative of a true decorrelation effect as 
predicted by Cronyn (1972) , although the decorrelation is not as severe as 
his model proposes. There is also no difference between the amount of 
decorrelation as a function of solar elongation. This is not unexpected, 
though, because our observations were made when the Crab nebula was be- 
tween 180° and 80° elongation and the source was never very near the sun 
where stronger effects are expected. 

Figures III-2 and III-3 are plots of the individual-run fluxes for 
111 and 74 MHz in 1973 February. They comprise our "visibility curves" 
for the compact source. Most of the scatter in the data points is attri- 
buted to random errors because the standard deviation is usually equal to 
the expected errors. However, we note that the large variations In flux 
at 111 MHz on 25 February (in figure III-2) are correlated on the long 
and short baselines so that the fringe visibility remains approximately 
constant. Therefore, some of the scatter in the average fluxes can be 
attributed to variations in the total flux from the pulsar on a time scale 
of minutes. On February 26 and 27 and for all of the 74-MHz observations 
(figure III-3) the pulsar was not "active" during the hour we observed it, 
and the scatter has the expected value. 

As illustrated in figure III-2 t values for the fringe visibility are 
surprisingly constant in spite of the rapid and 1 *rge variations in total 
flux. The unresolved and partially resolved fluxes followed approximately 
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Figure III-2. Total average flux of the Crab nebula pulsar. 111 MHz. Data 
is for three days In 1973 February. Each point is a three-minute average. 
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Figure III-3. Total average flux of the Crab nebula pulsar at 74 MHz. Data 
is for several days in 1972 December and 1973 February. Each point is a 
three-minute average. 
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the same trends on all times scales so as to keep the fringe visibility 
constant. Figure III-l shows this on a yearly time scale. Although some 
of the flux variation from month to month could be due to unknown system- 
atic errors, the fact that the long and short baselines usually vary in 
the same way argues for an interpretation of the changes which allows some 
intrinsic variability in the fulsar f s flux. The variability is discussed 
in section B-l of this chapter when the spectrum of the compact source is 
considered. 

3. Brightness distribution and angular size . The data presented in 
table III-l provide only a small range of points on the visibility curve 
because the sparse hour-angle coverage allowed little foreshortening of 
the long baselines. Interpretation of the visibility for any of the data 
can be done only by assuming a source brightness distribution. In the 
1971 March 144-MHz experiment the Crab nebula was observed at many dif- 
ferent hour angles so that each baseline was foreshortened, and a crude 
visibility curve is available. The data points are plotted in figure III-4. 
The fringe visibility scale was set by normalizing each observed correla- 
tion coefficient to the correlation observed on the shortest baseline. 

If the apparent angular size of the Crab nebula pulsar is due to in- 
terstellar scattering, then according to the theory we expect the size of 
the scattering disk to depend on the observing frequency to some power. 

For certain types of irregularity spectra, we also expert the source to 
have a Gaussian brightness distribution at all frequencies. Unfortunately, 
the 144-MHz data points are consistent with many brightness distributions. 
Superimposed on the data of figure III-4 are three visibility curves: a 

Gaussian (FWHM » 0.055 arcsec), a uniform disk brightness distribution 
(diameter 0.06 arcsec) and a point double (separation 0.04 arcsec). The 
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Figure III-4. Fringe visibilities of the Crab nebula pulsar observed in 
1971 March at 144 MHz. Model visibility curves for three brightness 
distributions are shown. 
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Gaussian curve fits as well as the other two curves, so we can only say 
that the data are consistent with a Gaussian brightness distribution. 

The observed fringe visibilities can be related to angular sizes by 
considering the Fourier transform relations. Assume the source's Gaussian 
brightness distribution is given by 

B(6) - exp [-4Hn2(e/e g ) 2 ] (III-l) 

where 0 = full width at half maximum (FWHM) of the Gaussian curve. The 

s 

two-dimensional Hankel transform of B(0) is the visibility function: 

y(b x ) * exp [-ir 2 b x 2 0 s 2 /4Jln2j (III-2) 

where b x * baseline length in wavelengths. We observed a value of y at 
some b, and we can calculate 0 from equation III-2 or read its value 

A 8 

from a graph of equation III-2 shown in figure III-5. The values of 0 

s 

obtained from the fringe visibilities in table III-l and figure III-4 are 
listed in table II1-2 and plotted in figure IXI-6, along with several other 
measurements for comparison. The errors in the determination of 0 g are 
derived from the measurement error in y by the formula for error propagation. 

Table III-2 

Apparent angular diameters of the Crab nebula pulsar 



Baseline ■ 

.aseline (MX) 

Y __ 

0 a (arcsec) 

196.5 

2550 

1.7 

0.74±0.06 

0.03±0.01 

144.3 

2700 

1.3 

0.65±0.15 

0.05±0.03 

111.5 

2550 

0.94 

0.72±0.03 

0.07±0.01 

74.0 

2550 

0.63 

0.32+0.02 

0.18*0.01 



Figure III— 5. Fringe visibility vs. angular size of source for a 
Gaussian model. Analytic forn for the curve is given in equation III-2 
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Figure III-6. Measured values of the apparent size cr 'he Crab nebula pulsar 
at several frequencies. Sizes ' eve been converted to Russian half -power 
diameter when necessary. The VLBI p tint at 26 MHz is by Mutel et al. (1976), 
the others are fro: thLs thesis. The references for the IPS points are: 

S Cronyn (1970a). Armstrong et al. (1973), 0 Bel’ & Hewi~h (1967). 

The line through the points has a slope of 2. 
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The fractional error is 


Ae s /0 s = <- 1 / An Y) * (A y/y) (HI-3) 

The value c. is much better determined if y is very small even if it 
has a large uncertainty. 

A weighted least-squares fit to a power-law wavelength dependence was 
made for the four values of 9 g in table III-2. An exponent of 2.0±0.4 
was ob - ' red , which is consistent with the wavelength-squared dependence 
found i. uauy scattering theories. The wavelength dependence of inter- 
stellar scattering is discussed in rare detail by Mutel et al. (1574). 

Tb°v have also collected many measurements of the apparent angular size 
of the compact source in the Crab nebula, and a fit to all of the data 
yields a power law oi exponent 2.05±0.25. 

B. Discussion 

1. Spectrum . The fluxes measured at our three frequencies give a 
spectral index of -2.9±0.4, which is the same value and error quoted by 
Rankin e£ al . (1970) from observations of the pulsar between 430 and 196 MHz. 
Our total flux measurement at 196 MHz agrees with the Arecibo flux measure- 
ment at that frequency, and our measured total fluxes at the lower fre- 
quencies therefore lie on the extrapolation of the spectrum of the pulsing 
flux. Figure III-7 give*, the soecttum of the compact source, including 
pu.aiir, IPS, and VLBI measurements. 

The usefu. .>.ess of a spectral index for this source can be questioned 
however, becaus® of the variability observed in both the pulsing and the 
total flux. Si jultaneouo measurements at several frequencies have shown 
there is no correlation between variations in the pulsing power at 430 and 
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Figure III-7 . Spectrum of Crab nebula compact source and pulsar. References 
for the spectral points are: $ Cronyn (1970a); ID Harris (1973); 0 Readhead 
and Hewish (1974); 0 Armstrong et al. (1973); 0 Erickson et al. (1972); 

O this thesis; 0 Rankin et al. (1970). 
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111 MHz (Rankin et al. 1974), which means that even the shape of the spectrum 
is variable* Changes in the pulsing flux of up to a factor of three in a 
5 w weeks time aa common at 430 MHz where all of the flux is pulsing, while 
at j.ower frequencies the changes in pulsing flux, though still apparent, are 
not so drastic. 

The behavior of the pulsing flux contrasts with our observations of 
the total flux. We observed up to 25Z variations in the total flux at 
111 MHz, and no significant variations (<15Z) at 74 MHz. This behavior 
also differs from that of other pulsars which seem to be highly variable 
at meter wavelengths (Huguenin et al. 1973). More monitoring ot the total 
flux is required to be able to derive statistical properties of the vari- 
ations which can be compared to the variations in pulsing flux. 

2. Small structure in the nebula . Several lines of evidence suggest 
that there could be small-scale emission structure (other than the pulsar) 
in the Crab nebula at low frequencies. The observations of Matveenko (1968) 
indicate there are two compact sources in the nebula. Using recent inter- 
planetary scintillation observations Armstrong et al. (1973) suggest that 
chere may be several small sources in the nebula. It is also possible that 
the wealth of detail present in the filaments, wisps, and continuum emis- 
sion at optical wavelengths might continue to low frequencies. High- 
resolution maps at centimeter wavelengths (Wilson 1972) reveal "ridges" 

-*hich are coincident with the filaments and have spectral indices approx- 
ini. tely equal to that of the continuum radiation. This information hints 
that these structures might also be detec ible at meter wavelengths. 

Although our coverage of the (u,v) plane is meager, we can draw some 
conclusions about small-scale emission structure in the nebula. A lack 
of fine structure is indicated by our flat visibility curves at 111 and 



61 


74 MHz (figures III- 2 and III-3). The fringe visibility as a function of 
hour angle at 111 MHz does not vary except for noise fluctuations. Also, 
there is no systematic trend in the time-delay and fringe-rate residuals 
which might indicate the presence of other sources. At 74 MHz the data 
spans two and one-half hours and again no changes in fringe visibility or 
residual parameters are evident. 

If there were fine structure besides the pulsar in the nebula at our 
frequencies, that is, if there were many small sources, the visibility 
curve should show large fluctuations with small changes in hour angle. 
Except for a fortuitous arrangement of sources we conclude that there is 
no st acture in the nebula on an angular scale of the order of the fringe 
spacing (0.3 arcsec at 74 MHz) or smaller. Similarly the smoothness of 
the short-baseline interferometer visibility curves also indicates a lack 
of structure with angular sizes smaller than 15 arcsec. More extended 
observations at these frequencies should be made in order to consider this 
problem in more detail. 

Another argument which supports the absence of structure in the nebula 
relies on data presented in the next chapter. Anticipating the results, we 
observe that all of the flux at 196 MHz is contained in the pulsing com- 
ponent, and that at 111 MHz the phase of the pulsar 9 s radiation agrees with 
that of the compact source. This is evidence that the pulsar is truly the 
only small object in the nebula, because if we were observing the coherent 
sum of radiation from the pulsar plus fine structure or other compact ob- 
jects, then part of the total flux would not be pulsing at 196 MHz, and the 
continuous radiation would have a different phase from that of the pulsar. 
This is not the case, and we conclude there are no compact objects or 
structure in the nebula other than the pulsar. 
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3. Discrepancy with IPS flux . Interplanetary scintillation (IPS) 
observations of the compact source in the Crab nebula at 74 MHz have been 
reported by Armstrong et al. (1973). They derive an angular diameter of 
0.18±0.Q1 arcsec at 1/e which corresponds to a half-power diameter of 0.15 and 
is consistent with our observations. The amount of scintillating flux 
which they quote is a factor of two higher than the total flux observed 
with VLBI techniques, although their value for the flux agrees with other 
IPS observations of the compact source by Harris (1973) and Readhead and 
Hewish (1974). The spectrum of figure III- 7 contains several points which 
illustrate the discrepancy between the IPS and VLBI fluxes. 

Armstrong et al. suggest there could be another source or sources in 
the nebula which could account for the extra flux they observe. The 
scintillating flux cot* Id come from any location within the telescope beam, 
and they would observe an incoherent average of the various source fluxes. 

Even though our hour angle coverage Is incomplete we have shown above that 
there is no "structure" or other sources In the nebula except for a pos- 
sible coincidental arrangement. 

If there Is only the pulsar which is small at 74 MHz they suggest we 
might be missing some of the flux due to decorrelation or other effects, 
but we believe we have accounted for all the flux. Decorrelation due to 
interplanetary phase fluctuations over the short 50-km baseline could not 
have occurred, especially since we Were observing when the Crab nebula 
was near the anti-solar direction. Ionospheric phase fluctuations were 
accounted for by broken coherence averaging when recessary. 

Accurate position measurements at low frequercies are impossible 
because of the unknown phase shifts imposed by tue ionosphere, but we 
never found a correlation peak at any position in the nebula other than 



63 


near the pulsar. The small residual time delays and fringe rates serve 
to confine the location of the compact source to the inner one-tenth of 
the nebula. 

The discrepancy between the IPS and VLBI techniques may have been 
resolved recently. Coles et^ al. (1974) obtain a good fit to the curve of 
scintillation index vs. solar elongation by using new theoretical calcula- 
tions and normalizing their data with the VLBI-meas ired flux. They con- 
clude that they now agree with our results* Further comparisons between 
VLBI and IPS results have been made by Resch (1974) using the data from 
these experiments and the IPS survey by Readhead and Hewish (1974). 

4 f Polarization. The instrumental polarization was very high, as 
shown in tables II-5 and II-6, and it was not well determined. Thus the 
treasured linear polarization of the pulsar's radiation has large uncertain- 
ties. Table I1I-3 lists the fraction of linearly polarized radiation for 
the total flux (m^,, from the short-baseline observations) and the partially 
resolved fluxes ( m , from the long baselines). The limits on the linear 
polarization at 111 and 74 MHz are due to the large error associated with 
the determination of the instrumental polarization. At 196 MHz, about 50% 
linear polarization was observed, which is consistent with the observations 
by Manchester et al . (1972). They observed a mean linear polarization of 
45% at 147 MHz and 55% at 281 MHz. Polarization of the pulsing power alone 
will be considered in the next chapter. 


Table III-3 


Polarization of total average flux 


Frequency 


ul 


”tr 


196 MHz 

111 MHz 
74 MHz 


0.40 


.40.10 

-0.30 


<0.15 

<0.15 


0.50 


40.10 


- 0.20 


<0.25 

<0,25 
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C. Summary of Chapter III 

Values for the total flux of the compact source in the Crab nebula at 
74, 111, and 196 MHz were presented in this chapter. The source was ob- 
served with both short- and long-baseline interferometers so that measure- 
ments of the fxi e visibility allowed angular diameters to be estimated 
at each frequency. The angular sizes scale as frequency-squared which is 
consistent with interstellar scattering theories. 

An absence of small structure other than the pulsar in the Crab nebula 
is indicated by the lack of structure in our observed visibility functions. 
The results of these measurements combined with those of the next chapter 
will be compared to predictions of interstellar scattering theories in 


Chapter V 



CHAPTER IV 


CRAB NEBULA PULSAR: PULSING FLUX RESULTS 

The results of our determinations of the Crab nebula pulsar* s pulse 
shapes are discussed in this chapter. The profiles were reconstructed 
from the Fourier components according to the method outlined in Chapter II, 
and further analysis is described in the first part of this chapter. Sample 
profiles are presented in the second and third parts of the chapter, and 
the implications of differences and similarities between the correlated 
and single-telescope pulse profiles are discussed. 

The single-telescope pulse shapes at several frequencies are shown 
in figure IV-1. which is taken from the paper by Rankin et al. (1970). 

The general pulse morphology 
changes with frequency because 
the components of the pulse 
which are distinct at high 
frequencies become blended 
due tc scattering. There are 
; components of the pulse, 
all of which contain com- 
parable amounts of energy. 

The precursor occurs just 

before the main pulse, and an interpulse is located approximately half 
way through the pulse. Between the main pulse and interpulse is a base- 
line component. 


73.8 MH* 



<30 MHt Nfr/Vtyy 

0 10 20 30 

MlLLtSCCQNOS 

Figure IV-1. Single-telescope 
pulse shapes. 
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A. General discussion of Fourier components 

1, Review of method . Besides the total time-averaged flux from the 

compact source in the Crat* nebula, we also considered separately the two 
components of the radiation: the pulsing (I p ) and non-pulsing or steady 

(1^) fluxes. As described in Chapter II the data tapes are processed with 
appropriate time-delay and fringe-rate offsets so as to make the interferom- 
eter respond to the nth Fourier harmonic of the pulse shape* The tapes are 
processed 2N+1 times to yield values of fringe amplitude and phase for n=0, 
±1, ±2, ... N, where N is the highest order harmonic which can be detected. 
The pulse shape is recovered by adding the complex components, including 
n^O, in a Fourier sum, as discussed in detail in Chapter II (pages 22-27), 

2, Number of components . For a simple single pulse shape (ignoring 
the Crab nebula pulsar f s interpulse for now) the amplitude of each Fourier 
component decreases monotonically with increasing harmonic number, but the 
rate of decrease depends on the pulse width. A sharp pulse (such as the 
196-MHz pulse) has a broad Fourier transform and high-order components 
have appreciable amplitude, whereas a broad pulse (such as at 111 MHz) 
has a narrow Fpurier transform so that the amplitudes of the components 
fall rapidly with higher harmonic number. With comparable sensitivity at 
the two frequencies, we therefore expect to be able to detec - more com- 
ponents at 196 MHz than at 111 MHz. In practice, we usually detected five 
Fourier components in one run (one three-minute average) at 196 MHz. Since 
the pulses do not have simple shapes, at 196 MHz the amplitude of the n“2 
component is higher than that of n*l because of the strong influence of 
the interpulse. At 111 MHz three Fourier components were usually detect- 
able in one run, although the actual number varied with changes in the 
pulse shape and level of pulsing power. 
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As discussed in Chapter I, we had expected to observe a narrower 
pulse with the long-baseline interferometer than with the short baseline 
because of the scattering geometry. This means we should have detected 
more Fourier components in the higher resolution data. For example, we 
might have observed a pulse similar in shape to the 196-MRz pulse when 
we observed with the long~baseline interferometer at 111 MHz. This turned 
out not to be the case, as the pulse shapes on all baselines are similar. 
This fact and its implications are discussed in later sections. 

Because the pulse waveform is a real quantity, we expect the +n and 
-n components to have equal amplitudes and opposite phases. Due to the 
presence of noise, this is not always the case. Corresponding components 
hardly ever have equal amplitudes, and the phase noise is typically about 
one-half radian. Occasionally, for high values of n, one or the other of 
a component pair will not be detected. On some tapes, perhaps, components 
up to n=8 at 196 MHz could be detected, but to prevent waste of computer 
time we processed the tapes for only five components. Similarly, experi- 
ence showed that n=3 was the highest component number to be used at 111 MHz 
because higher order components were rarely detectable. 

Tests were made in which ** e component amplitudes and phases were 
perturbed by various amounts in o^der to investigate the effect on the 
profile. It was found empirically that the phase could be perturbed by 
up to one-half radian and che amplitudes coaid be adjusted by about 50% 
without distorting the original profiles beyond recognition. Correct 
amplitudes but u rrresponding random phases do not result in a "believable 11 
profile. The waveforms at 111 MHz observed at Arecibo under the Crab 
nebula pulsar surveillanr j program are analyzed into six Fourier harmon- 
ics (Counselman and Rankin 1 72) . twice as many as we could detect. This 
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is presumably due to the fact that the components with n > 3 are only 
marginally detectable. Since these components are very weak, their 
presence or absence does not r;ffect the pulse profile significantly. 

3. Lack of components at 74 MHz . At 74 MHz we never detected any 
Fourier components with the long- or short-baseline interferometers, al- 
though the single-telescope waveforms observed at Arecibo can usually be 
analyzed into at least one Fourier component. This discrepancy is not 
readily understandable; it could possibly be due to several effects which 
a^e co' 'idered next. 

One difference between the two observing methods is the relative 
sensitivities of the Mark I system and the Arecibo receivers. Table IV-1 
gives the parameters for the two systems; the Arecibo parameters were 

Table IV-1 

System parameters for Crab nebula pulsar observations 


System 

Freq. 

(MHz) 

Bandwidth 

(kHz) 

Integr . time 
(min.) 

VLB I Mark I 

all 

350 

3 

Arecibo 

196 

4-10 

18 

Arecibo 

111 

2-10 

18 

Arecibo 

74 

1-5 

18 

taken from the paper by Rankin j 

et al. (1970). 

These and other system 

parameters from cables 

II-2 and 

II-4 were used 

to calculate the sensi- 

tivity, when observing 

the Crab 

nebula, for „< 

ecibo (A) and for each 


interferometer (double letters). Table IV-2 presents these numbers. 
The zero-length baseline (AA) is formed by "auto-correlating" and de- 
dispersing the Arecibo Mark I data tape. This procedure was discussed 
in Chapter II, section B-3 (page 27). 
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Table IV-2 

System sensitivities in Jy 
for observations of the Crab nebula 


Frequency 


(MHz) 

A 

AA 

AN 

AS 

NS 

196 

2.8 

1.7 

1.7 

2.3 

2.3 

111 

4.1 

2.0 

4.2 

3.0 

3.4 

74 

5.1 

2.4 

3.6 

4.5 

7.0 


Our short-baseline sensitivity is nearly the same as that of Arecibo 
but no Fourier components were detected on a few test runs. This lack 
of detection is not due to the flux being below the detection limit since 
on these rims the Arecibo system recorded between 15 and 25 Jy of pulsing 
flux (Rankin et al. 1974) . With th' 1 long baselines , the interferometer 
partially resolved the total flu * and the amplitudes of the Fourier com- 
ponents are reduced by this factor also because the effect of a sharper 
pulse with high resolution does not occur. This puts even the n»l com- 
ponent well below the detection limit on the long baseline. 

Table IV-2 shows tint the zero-baseline interferometer formed using 
the Arecibo data is the most sensitive; this is due to the relative* ’* 
Mark I bandwidth. We did detect the first Four .er component j an 
amplitude of about 3±1 Jy on one data tape from 1973 February 28. On that 
day the Arecibo system detected 15±6 Jy of pulsing flux (Rankin e£ al . 

1974). Even considering the large errors, our marginal detection is not 
entirely consistent with the Arecibo observations. 

One possible problem Involves curvature of the pu'<. ^s •» ■ >ss our 
passband since, as shown in Chapter II, we approximate the dispersion of 
the pulses as a straight line. The effect of the curved pulse was estimated 
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by calculating the difference between the pulse time-of-ar rival at the 
top and bottom of the passband using both the straight-line approximation 
and the true dispersion relation. At 74 MHz the difference between these 
two times is about 3 msec. This corresponds to about one-half radian of 
phase difference when we integrate across the passband, resulting in a 
signal degradation by a factor about 0.7. This would decrease the zero- 
baseline interferometer’s sensitivity but only to about that of the Arecibo 
system. 

Several runs were processed in searches for jurier components , but 
no general production runs were made because of the low signal le/el. In 
order to draw better conclusions about possible difficulties with the de- 
dispersing method at very low frequencies more runs should be processed to 
confirm the lack of detection; this processing ic in progress, 

4. Epoch determination . The phase of each Fourier mponent is 
referenced to the time at which tne offset oscillator used in the pro- 
cessing was started. r Tiis time was tha same for all values of n. This 
means that the peak of the reconstructed pulses occurs at a phase which 
can be referred to UT. We have compared pulse arrival times obtained 
this way to those measured at Arecibo and the epochs agree to within a 
few percent of the pulse period. Pulse shapes an be coherently aver- 
aged using the epoch information if accurate times-of-arrival can be 
predicted. 

B, Results at 196 MHz 

1. Sample profiles and Fourier components . The five Fourier com- 
ponent amplitudes and phases for each three-minute run at 196 MHz were 
added in Fourier sums, Interpolated, and the resulting smooth profile was 
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plotted. Figure lV-2(a'> shows typical pulses from the short and long 
baselines and Arec’bo. The component amplitude.. .aid phases ior the pro- 
files are in the next two diagrams: Figure IV-2(b) is a plot of the 

Fourier component amplitudes, normalized to the n=0 component amplitude, 
against harmonic number; figure 17-?. (c) is .. plot of the phases, 

relative to the phase of n*0, for the shor»--ba. ; .elir.e profile. The phases 
of each component pair (+n and -n) aie symmetrical .'bout the phase of n*0 
with an rms deviation of 9?4 for this r^ Appendix E contains a catalogue 
of 196-MHz average p -.lse shapes and a listing of the Fourier components and 
other relevant quantities. 

2. Discussic n. Both of the interferometer profiles and the Arecibo 
profile are similar and the amount of resolution Is slight. All of the 
flux is pulsing (I,p = Ip), and the baseline level set oy the interferom- 
eter is the same as that of the tingle-telescope observation. The reason 
for not using 196 MHz for the second set of experiments was that the in- 
terferometer provided no unique information compared to a single telescope. 

Since all of the power is contained in the pulsing component, the 
linear polarization of the totrl flux xs the same as that of the pulsing 
flux. The traction of linear poi irizotion, m, is —50%, ar< in Table III-3. 

We have not attempted to determine m separately for the main pulse, and 
interpulse because of th*’ high instrumental polarization. 

C. Res ul ts at 111 MHz 

1. Sample profiles and variation of Ip/X^, At 111 MHz tt* pulsing 
(Ip) and steady (Ig) components of the pulse shape contribute about equally 
to the total flux from the compact source (I , - I r — I_/2). Thus the ratio 
Ip/I T is approximately equal to one-half, although it was iound to vary with 
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FLUX (Jy) 




HARMONIC NUMBER 



Figure IV-2. (a) Typical pulse profiles, 196 MHz. (b) Amplitudes, 

relative to n«*0, of Fourier components of the pulses in (a), (c) Phases, 
relative to n»0, of Fourier components of the short-baseline pulse in (a). 
Vector length indicates component strength. 
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time. In the data from the first five experiments, Fourier compon nts 
up to n=3 were usually detected. The resulting pulse shapes are 01 .4*- 
than the 19 6 -MHz pulse shapes and they also have an offset baseline bias 
and continuous non-pulsing power. Pulses from a typical run in 1971 No- 
vember are shown in figure IV-3(a). The comparison Arecibo pulse has been 
raised above the baseline to coincide with the short-baseline interferom- 
eter pulse. Figure IV- 3(b) is a plot of Fourier component amplitude vs. 
harmonic number; figure IV-3(c) is a polar plot of the phases. The aver- 
age (signed) deviation of the phases from symmetry around n=0 is 9?8 and 
the rms deviation is 24“ 6. 

.The value of Ig is taken as the minimum of the pulse profile. The 
total flux Xj is the amplitude of the n=0 component (A q ) , and Ip is calculated 
as Ip = 1^. - Ig. This underestimates Ig because some "ringing” of the 
Fourier harmonics is expected and the minimum of the profile may be below 
the true base level. This is evident on the 196-MHz profiles where the 
amplitude is sometimes negative. For the first set of five experiments 
(1971-2) the average value of Ip/Ip was 0.6 for both long- and short-base- 
line pulse profiles. The ratio of Ip/ Ip is 0.72 for the short-baseline 
profile of figure IV-3. 

For the last two experiments most of the data revealed only two 
Fourier components, and Ip/Ip ~ 0.4 for all the profiles. Figure IV-4(a) 
is a reconstructed pulse observed in February 1973. The pulse shape is 
clearly broader than that in figure IV-3. Figures IV-4(b) and IV-4(c) 
give the component amplitudes and phases relative to n*0. Appendix F is 
a catalogue of 111-MHz average pulse shapes from all experiments, and a 
list of the components and values of Ip/ Ip. 
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Figure IV-3. (a) Typical pulse profiles, 111 MHz, 1971 November. 

(b) Amplitudes, relative to n*0, of Fourier components of the pulses 
in (a), (c) Phases, relative to n-0, of Fourier components of the 

short-baseline pulse in (a). Vector length indicates component strength. 
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Figure IV-4. (a) Typical pulse profiles. 111 MHz, 1973 February. 

(b) Amplitudes, relative to n«0, of Fourier components of the pulses 
in (a), (c) Phases, relative to n M 0, of Fourier components of the 

short-baseline pulse in (a) . Vector length indicates component 
strength. 
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Figure IV-5 is a plot of the average value of Ip for each observing 
day, and also the average value of Ip/I^. Table IV-3 lists the values 
in the plot and the fringe visibility of the pulsing flux, y • The ex- 
pected uncertainty in Ip is quite large, about ±15%, because Ip is de- 
rived from a sum of 2N+1 amplitudes and phases which each have an error. 

The scatter about the mean value of Ip is about ±20%, only slightly higher 
than the predicted value. 

The total amount of pulsing power was determined in a variety of 
ways: with the short-baseline cross-correlation interferometer, by auto- 

correlating the Arecibo tape to form a zero-baseline interferometer, and 
from the single-telescope observations made at Arecibo. The results of 
each method are plotted in figure IV-5. All these methods give approxi- 
mately the same value of Ip, although the short-baseline correlated flux 
is systematically lowest by about 20%. This is somewhat disturbing since 
it occurs only for the cross-correlation observations. However the large 
error bars allow the several flux values to be consistent. The Mark I 
auto-correlation data agrees well with the Arecibo fluxes which were ob- 
served simultaneously with the VLBI experiment. 

2. Polarization . The polarization of the total flux at 111 MHz 
was presented in table III-3. With the pulse profiles in hand we can 
separate the pulsing and non-pulsing flux and determine the polarization 
of each component separately. Table IV-4 lists nip (the fraction of linear 
polarization in the pulsing power) and m^ (fraction in the total power, 
from table III-3). The values for nip are consistent with the results of 
Manchester et al . (1972) who observed nip * 35% at 114 MHz. There is a 
formal difference between the value of the polarization of the pulsing 
and total J luxes, but within experimental uncertainty the two polarizations 


are the same 
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Figure IV-5. Average values of pulsing flux (I. 
and the ratio of pulsing to total flux (I_/I T ). 




Pulsing flux of the Crab nebula pulsar at 111 MHz 


78 






m 




CO 

o 

m 


rH 

o 


o 




rH 

*— l 

o 


• 

• 


• 




• 

• 

• 


O 

o 

I 

o 

1 

i 

1 

o 

o 

o 

o. 

+1 

+i 

\ 

+i 

1 

i 

1 

+1 

+1 

+i 

>- 

On 

o 






CM 

rH 

m 


r*. 

00 


m 




00 

00 



• 

• 


• 




• 

« 

• 


O 

o 


o 




o 

o 

o 

• 1 


*■3- 

ON 

CM 

<r 

i 

1 

in 

NO 

CM 

O 


iH 


rH 




rH 

rH 

CM 

35 1 












O 

vO 

r^. 

rv 

ON 



o 

rH 

O 


• 

• 

• 

• 

• 



• 

• 

• 

Pd 

<r 

m 



CO 

i 

1 

in 

>3” 

CO 

PU 

+i 


+1 

+1 

+\ 

i 

1 

+i 

+1 

+1 

M 

m 

ON 


O 

00 



rH 

ON 

CO 


♦ 

• 

• 

• 

• 



• 

• 

• 


CM 

00 


rv 

CM 



CO 

CO 

m 


CM 

rH 

i— f 

rH 

H 



rH 

rH 

rH 


ON 

NO 


co 

VO 




CO 

ON 

Pd 

o 

i— 1 

o 

iH 

O 



iH 

rH 

O 

H 

• 

• 

• 

• 

• 



• 

• 

• 

M 

O 

O 

o 

O 

o 

1 

1 

o 

o 

o 


•fl 

+1 

+i 

+1 

+1 

1 

1 

+1 

+1 

+1 

Pd 

NO 

CO 

iH 

ON 




00 

00 

rH 

PH 

VO 

NO 

NO 

m 

m 



sr 

<r 

m 

►H 

• 

• 

• 

• 

« 



• 

• 

• 


o 

o 

O 

o 

o 



o 

o 

o 

• 

m 

CO 

1 

<r 


rH 

c0 

in 

•H 

rH 

O 






rH 



rH 

rH 

z 












co 

NO 




NO 

O 

<r 

CM 

CO 


• 

• 


• 


• 

• 

• 

* 

• 


NO 

rH 

1 

CO 


CO 


-3“ 

m 

CO 

Ph 

+1 

+1 

1 

+1 

| 

+1 

+i 

HH 

+i 

+1 

M 

vO 





ON 


O 

rH 

<r 


• 

• 


• 


» 

• 

• 

• 

• 


00 

co 


On 


rH 

co 

NO 


o 


CM 

CM 


CM 


rH 

rH 

rH 

rH 

CM 


CO 



On 


rH 

O 

CM 

rH 




o 


O 


rH 

rH 

rH 

rH 

O 

H 

• 

* 


• 


• 

• 

• 

• 

• 

M 

o 

o 

1 

o 

l 

O 

o 

O 

o 

o 


+1 

+1 

1 

+1 


+1 

+1 

+1 

+1 

+t 

PM 

ON 

co 


rH 



o 

On 

CM 

NO 

H 

uo 

m 


NO 


co 

>3“ 

co 


<r 


♦ 

• 


• 


• 

• 

• 

• 

• 


o 

o 


o 


o 

O 

O 

o 

o 


<3* 

o 

Mf 

00 

m 






CM 

CM 

CM 

CM 

CM 






A 

ON 

1 

CO 

1 

CM 



in 

NO 


CM 

rH 

CM 

CM 

CM 

00 

ON 

CM 

CM 

CM 

• 

> 

• 

u 

• 

C 

• 

X* 

• 

u 

• 

O 

• 

o 

• 

rO 

• 

X* 

• 

x> 

O 

<u 

3 

0) 

cd 

0) 

4) 

0) 

0) 

Q> 

55 

P4 

*o 

fe 

X 

Q 

o 

pH 

fe 

PH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CO 

CO 

CO 


r* 

Cs 








ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

rH 

rH 

H 

H 

rH 

rH 

rH 

rH 

rH 

rH 



79 


Table IV-4 

Polarization of total 
and pulsing flux at 111 MHz 


Fraction 

“P 

“T 


Total flux 


Resolved flux 


0.27 


,+0.10 


- 0.20 


0.22 


,40.15 

- 0.20 


<0.15 


<0.25 


3. Dependence of visibility on pulse phase . The predicted de- 
pendence of fringe visibility on pulse phase, illustrated in figure 1-3, 
was not observed. The pulse shapes observed on both long and short base- 
lines appear to be similar, in that the pulse derived from the partially 
resolved flux is uniformly reduced in amplitude but not changed in shape. 
Values of Ip/I^ were calculated separately for the pulses observed with 
long and short baselines and, to within the rather large error limits, 
are equal. Also, figures IV- 2b, 3b, and 4b indicate there is no system- 
atic trend toward having y increase (approach unity) for higher Fourier 
components . 

A quantitative test was made by computing the mean values for the 
flux of each Fourier component using all the available data. (The In- 
dividual components are tabulated in Appendices E and F.) Then the fringe 
visibility was calculated for each component at both 111 and 196 MHz, and 
the results are plotted in figure IV-6. The error bars represent the 
standard deviation of y as calculated from the standard deviations of 
the mean fluxes. The figures show that each component has approximately 
the same visibility as the n*0 component, and within the errors there is 
no trend toward unity. All of the phenomena mentioned in this section 
are of course inter-related, and they all give useful ways of indicating 
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Figure IV-6. Figure visibility vs. harmonic 
number for (a) 111 MHz and (b) 196 MHz. 

the lack of strong evidence for a source which expands with pulse 
phase. 

4. Identification of pulsar with compact source . The arguments 
which support the identification of the low-frequency compact source with 
the pulsar in the Crab nebula were sunmarized by Erickson e_t al . (1972) 
and also by Lang (1971b). High-resolution VLBI observations of the Crab 
nebula can be used to distinguish between the compact source and the pul- 
sar if the phases of the radiation are known. The correlated flux obtained 
through processing the data as for a continuum source is due to the compact 
source, and the phase of these fringes provides the position of the com- 
pact source. The fringe amplitudes and phases obtained for n ^ 0 obviously 
are due to the pulsar, and the average phase of the +n and -n components 
will give the position of the pulsar. Only if the compact source is in 
the same position as the pulsar will the fringe phases be the same. The 
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phases of the Fourier components which were detected at 121 MHz were used 
by Erickson et al. (1972) to identify the pulsar with the compact source. 
The data presented in this thesis extend their results and tighten the 
limits on the identification. 

The method they used, and which is used here, is the measurement of 
differential phase, which is a powerful tool for determining accurate 
relative source positions. If the sources on which the differential mea- 
surement is to be made are in the beam of the telescopes at the same time 
then all instrumental phase deviations, including differences between lines 
of sight through the atmosphere and ionsophere, are greatly reduced. The 
difference in fringe phase between the sources is then a direct measure of 
their position difference. The test is very sensitive because we can mea- 
sure a difference in phase to a fraction of a fringe. 

The problem is formulated as a null test: the signed deviations of 

the average phase of ea : : pair (+n and -n) of components from the phase 
of the corresponding n=0 component are summed and averaged to yield 
If - 0, then the pulsar would occupy the same position in tlje sky as 

the compact source. The values of are listed in the tables of Ap- 

pendices E and F, for individual runs, and were averaged fo* the values 
listed in table IV-5. The averages were formed for both frequencies. The 

Table IV-5 


Average phase differences between pulsar and compact source 



196 MHz 2.514.0 0.811.3 41 205 


111 MHz 


-1.413.6 


0.812.0 


98 


241 
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values were computed using all the long-baseline data, quoted 

error is the standard deviation of the mean. At 111 the ? cage phase 
difference implies that the pulsar and compact source a^e coincident to 
about 0,8 milli-arcsec in a position angle along the baseline. This cor- 
responds to a linear distance limit of —2 a.u., assuming 2 kpc for the 
distance to the Crab nebula (Trimble 1973). 

As a test of the method, the same procedure was followed using the 
196-MHz Fourier components. At this frequency we know that all of the 
power is pulsing because the baseline level of the pulse prof Lies is at 
zero, so must be zero. The error gives an idea of the noise level. 

As the numbers in table IV-5 show, the result and error are approximately 
the same as at 111 MHz. These results were used in Chapter III, part B-2 
as evidence against small structure in the Crab nebula other than the 
pulsar. 

D. Summary of Chapter IV 

The body of data presented in this chapter consists of the Fourier 
components of the Crab nebula pulsar's radiation at 196 and 111 MHz. 
Analysis of the amplitudes, phases, and pulse shapes has shown that the 
shapa changes with time and that the baseline offset contains a variable 
fraction of the total flux. We have seen that the fringe visibility of 
all components of the radiation is not a function of pulse phase. Finally 
the pulsar has been undoubtedly identified with the compact source, and is 
coincident with it to within 2 a.u. at the distance of the nebula. The 
next chapter will use these results and those of the previous chapter to 
examine various models and theories. 



CHAPTER V 


MODELS FOR THE CRAB NEBULA AND INTERSTELLAR MEDIUM 
A, Review of observed phenomena 

Many of the spatial, temporal, and spectral features which we observe 
in the radiation from the Crab nebula pulsar are due to the influence of 
the media through which the radiation passed. Models for the radiating 
mechanism, and of the neular, interstellar, and interplanetary media 
should account for all the observed phenomena, a brief review of which 
follows. 

1. Spatial features . As described in the previous two chapters, 
the spatial features are uniquely observed via VLBI and T " " This thesis 
presents consistent observations of apparent angular «i. . a function 
of frequency. The brightness distribution has been observed at one fre- 
quency (144 MHz) and is consistent wit.* a Gaussian. The fringe visibil- 
ities which we measured have been interpreted in teims of Gaussian bright- 
ness distributions and yield the angular sizes listed in table III-2 
(page 55). The fringe visibilities measured for the total flux apply 
to all the components of the radiation: there is no systematic dependence 

of fringe visibility nn pulse phase, and the pulsing flux has the same 
visibility as the total flux. The technique of interplanetary scintil- 
lations can also be used to study angular structure of radio sources. 

Early IPS observations at 81 MHz (Bell end Hewish 1967) were used to 
estimate a size which is consistent vith our observations and recent 
IPS observations at 74 MHz (Armstrong elt ai. 1973) derive the same size 
as was measured with our VLBI experiments. 
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2. Temporal features . Many aspects of the temporal features of 
the radiation are being monitored with the twice-weekly observing program 
at Arecibo. During time intervals from weeks to months every observable 
parameter appears - vary in a non-systematic way. The pulce shape, dis- 
persion measure, amount of pulsing power, pulse period, and the period 
time-derivative all vary (Rankin and Counselman 1973). In the seven VLBI 
experiments, we observed in addition that the total flux and the ratio of 
pulsing to total flux are also variable. No temporal changes were obser- 
ved in the angular structure of the radiation. 

3. Spectral features . Narrow-band spectral features have not been 
observed for the Crab nebula pulsar, probably because the scintillation 
bandwidth of the features, predicted from the pulse broadening (~15 msec), 
is only about 10 Hz at 111 MHz. If a filter this narrow could be made it 
would be difficult to achieve an adequate signal- to-noise ratio. For 
jther pulsars, Rickett (1970) observed a relation between the dispersion 
measure and the maximum bandwidth over which scintillations are correlated. 
This implies ar. association between the "dispersing" electrons a-d the 
"scattering" electrons. They may be the same electrons, meaning that 
electron density irregularities occur wherever there are electrons. This 

is the basis for the belief that the scattering is truly an interstellar 
phenomenon and that it occurs all along the path. 

Several models have been proposed which attempt to explain these ob- 
served spatial, temporal, and spectral features, but so far none has suc- 
ceeded in predict! lg all of the observations. There are one or more as- 
pects of each theory which, considering especially the new data presented 
in this thesis, do not agree with the observations. In sections B through 
E which follow, we discuss four models and point out their virtues and 
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deficiencies. In section F, we have combined the best features of each 
model into yet another model which does in fact explain the observations. 
More refined observations in the future will require alterations to these 
theories . 

B. Thin screen model 

The first attempt to explain what happens in the interstellar medium 
to affect pulsars* radiation was the "thin screen model." This model was 
briefly discussed in the introduction and now will be considered in more 
detail. The discussion draws upon the work of Lang v.l971a), Cronyn (1970a 
and b) , and Lovelace (1970). 

1. Assumptions , (a) The artifice of the thin screen can be used 
if it is assumed that the density irregularities in the medium are iso- 
tropic and homogeneous; in view of Rickett's work this may be a fair 
assumption. The irregularities are collapsed onto a plane approximately 
midway between observer and source; the screen in thin enough that it 
produces only phase deviations in the wavefront while the amplitude remains 
unchanged, (b) Another assumption is that geometrical optics is a valid 
approximation so that ray tracing can be used, (c) Finally the time de- 
lays caused by variations in the group velocity through the density ir- 
regularities are ignored because they are much smaller than the delays 
due to geometrical effects. This last assumption is valid only for 
Gaussian and shallow power-law spectra of irregularity sizes, however. 
Cronyn (1973) specifically discards this assumption and proposes that 
the group delays are aa large as geometrical delays if a steep power-law 
spectrum of sizes is used. This important model will be discussed later. 
Jokipil (1973) notes that the thin screen is really only a convenient 
contrivance and not the limit of a more rigorous theory. 
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where q = L/ (Z - L) , (V-2) 

Z = total path length from source to observer, 

L = distance of observer from the screen. 

The quantity q is a proportionality factor which Indicates the location 
of the screen. It is the ratio of the distances of observer and source 
from the screen, and if the screen is halfway between observer and source 
q * 1. Note that as the screen is moved closer to the source q becomes 
very large, thus allowing a large time delay to be associated with a very 
small angular size simply due to the geometrical configuration. 
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The rms scattering angle is 0 q , defined through the two-dimensional 
brightness distribution: 

B(8) = exp - (6 2 / 26 o 2 ) . (V-3) 

6 is related to measured 0 by 
o s 

6 = 2(2*n2) 1/2 6 = 2.3 0 . (V-4) * 

S O o 

The time delay is the e-folding time of the exponential decay of 

the impulse response of the thin screen. The impulse response is also 

called the "probability density function" (pdf) for the scattering time 

—x 

delay t. For a single thin screen, the pdf has the form e which is 
plotted in figure V-2. 

The convolution of the 
impulse response with the 
emitted pulse shape yields the 
predicted pulse shape. For the 
calculations which were made 
here, the emitted pulse was 
represented as two delta functions: one for the combined main pulse, 
precursor, and baseline component; and one for the interpulse. The 
second delta function followed the first by 13 msec and had one-third 
of the total energy. The period was taken as 33 msec. This model is 
taken from Cronyn (1970a). A different model in which the delta func- 
tions had equal energies and were separated by 10 msec (combining the 
main pulse with the precursor and the interpulse with the baseline com- 
ponent) gave the same results to within a few percent for the calculation 
o f the fraction of pulsing flux. A more complicated model which would 



Figure V-2. Impulse response of 
thin screen. 
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include all four components separately was not tried, but for the purpose 
of calculating Ip/I^ such a model would probably not change the results 
significantly. Figure V-3 is the predicted pulse shape assuming ~ 15 
msec at 11Z MHz. 

The convolution of the pdf 
with a train of pulses yields 
the "steady state" observed 
power as a function of time, 

tlflUK) 

including any flux which is Figure V-3. Predicted pulse shape 

for thin screen. 

delayed enough to appear as 

unpuised radiation. Evaluating this function at the pulse period (the 

minimum of the curve) gives the amount of flux which is in the steady 

source, thus accounting for the contributions from the exponential tails 

of all previous pulses. For the thin screen ->del f figure V-4 is a plot 

1.0 

of the fraction -of the total 
flux which is still contrined ** 

in the pulsing component, as <u 

►- 
> 

a function of the character- a4 

istic tine delay. At a time 

0.2 

delay approximately equal to 

do 

the pulse period only about 

Figure V-4. Predicted I p /I T for 

one-third of the flux is still thin screen model. * T 

pulsing. 

3. Comparison to observations . The pulse shape of figure V-3 is 
obviously not in good agreement with the observed pulse shape because of 
the discontinuities. The exponential decay in the latter part of the 
pulse is observed however. 
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We observed Ip/I^ ~ 0.6 in 1971-2; from the graph of figure V-4 this 
corresponds to a time delay of 15±3 msec. If we let q = 1 in equation V-l 
and take the distance to the Crab nebula (Z) as 2 kpc (Trimble 1973) we 
calculate 0 g = 0.18±0.02 arcsec. This is more than a factor of two larger 
than the size which we observed, 0.07±0.01 arcsec. Alternatively, the 
observed angular size implies a time delay of only 2.3 msec. 

The calculated 0 g could be changed by making q ~ 6 and thus placing 
the screen at a distance L ~ 1.7 kpc from the earth. There have been sug- 
gestions that scattering may occur primarily in spiral arms which would 
constitute single screens. There is no strong spiral feature at 1.7 kp*_ 
in the anticenter direction, although the structure in this part of the 
Galaxy is not well known. 

Before our measurement of 0 g at 111 MHz, the size most often quoted 
was that derived from 81-MHz IPS observations by Bell and Hewish (1967) 
of 0 g = 0.2±0.1 arcsec which scales as wavelength-squared to 0.11±0.05 arc- 
sec at 111 MHz. If we assume this is a diameter at 1/e and convert to a 
half— power diameter, we obtain 0.09±0.05 which is consistent with 
our observations but, again, not with the observed pulse broadening. 

More recently Readhead and Hewish (1974) published an angular diameter of 
0.2510.1 arcsec, which is close to the previous value. The size derived 
from other IPS observations by Armstrong et al . (1973) at 74 MHz also 
agrees with our observations. 

Another problem with the thin screen treatment is that, to within 
our errors, we observed no change of visibility across the pulse. There 
is no way to avoid this phenomenon with the thin screen model since it is 
due to geometrical effects and derives from the assumptions mentioned 
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above. More accurate observations of correlated pulse shapes should be 
made before this effect can be ruled out completely, however. 

Besides the discrepancy between the temporal and angular broadening, 
a major problem with the thin screen theory is the concept: it is too 

simple. We know that between the earth and the Crab nebula there are 
many parsecs of interstellar matter, and it seems unlikely that the ef- 
fects along the entire path can be represented by a single thin screen. 

In fact the basic assumption of the thin screen approximation, that the 
medium is isotropic and homogeneous, is undoubtedly false. However, this 
approach was useful for the Crab nebula pulsar as a first step toward 
explaining the observations. 

For other pulsars the thin screen approximation also appeared to 
work well. Lang (1971c) showed that fits to the pulse shape predicted 

with the thin screen theory are generally quite good. Using more sen- 

sitive data on these same pulsars, though Co unselman (1974) showed that 
the residuals from a fit to the predicted pulse shapes for two pulsars 

are small but significant* Now that more accurate and extensive obser- 

vations of pulse shapes and angular sizes are available a more sophis- 
ticated theory than the single-screen model is warranted. 

C. Two-screen model 

1, The mod al. The next step toward improving agreement between ob- 
servations and theory was made by Counselman and Rankin (1972) who pro- 
posed that two scattering screens affect the radiation of the Crab nebula 
pulsar: the Crab nebula Itself and the interstellar medium. In a later 

paper, Rankin and Counselman (1973; hereafter called R&C) suggested that 
the nebular screen is variable in that sometimes it would be a stronger 
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scatterer than the interstellar medium, other times its effects would 
disappear, and often it would be comparable in scattering strength to 
the interstellar medium. R&C point out that the average electron den- 
sity and number of density irregularities in the 2-kpc path between the 
earth and the nebula should be constant in time, and they attribute the 
observed temporal variations in dispersion and pulse shape to changes 
taking place within the nebula. 

The identification of the nebula as a scattering screen was made 
because of other activity associated with the nebula and pulsar. About 
the time of the first observed large increase in dispersion measure (1969 
September), a major timing irregularity ( n glitch ,f ) also occurred. This 
activity was also accompanied by movement of the optically-observed wisps 
near the pulsar (Scargle and Harlan 1970). No such correspondence of 
these phenonema has occurred again, and it is likely that the 1969 Sep- 
tember events were chance coincidences. Nevertheless the concept of 
the nebula as a scattering screen is attractive. The nebula surrounding 
the pulsar must have some unique effect on the pulsar’s radiation, and 
the influence of the nebula should be assessed even if it can be shown 
to be negligible. This question is considered in section F. 

2. Theory . Functions similar to those plotted for the single- 
screen model will be developed now so that we can compare our observa- 
tions to the predictions of the two-screen model. The model of R&C is 
shown in figure V-5. They assigned 6.2 msec of time delay (at 100 MHz) 
to each screen, thus choosing 
an epoch when the nebula and 
interstellar medium were of 
comparable strength. The 
scattering pdf for two thin 
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screens with equal time delays (t^ ■ 

"“X 

a function of the form xe , 
which results from the con- 
volution of two exponential 
functions. 

The pulse shape result- 
ing from the convolution of 
this pdf with the model emit- 
ted pulse is shown in figure V-7, for 
t q - 6.2 msec. This predicted 
pulse shape is in much better 
agreement with observations 
than the shape predicted using 
a single thin screen (figure 
V-3) • 

The relative amount of 
pulsing power as a function of is 
equal time delays was used 
for simplicity. Comparison 
to figure V-4 shows that the 
pulsating component is cut 
off more efficiently with 
two screens than with one 
screen. Note that each of 
the two screens in the two- 
screen model has a time delay 


) is shown in figure V-6; it is 
o 
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Figure V-6. Impulse response 
of two-screen model. 


a frequency near 100 MHz with 
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Figure V-7. Predicted pulse 
shape for two-screen model. 


plotted in figure V-8. The pdf for 



Figure V-8. Predicted Ip/I_ for 
two-screen model. 


so that the scales for comparing figures V-4 and V-8 are different by 
approximately a factor of two. 
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If each screen has a Gaussian irregularity spectrum there should also 
be a decrease of visibility across the pulse as in the single-screen case, 
although the effect is not as pronounced here because the pulse is already 
partially broadened when it reaches the interstellar screen. The visibil- 
ity function for pulsar radiation scattered by two thin screens is dis- 
cussed in more detail by Cronyn (1974) . 

3. Comparison to observations . During all of our observations the 
time delays of the two screens were almost certainly never equal. The 
time delay which is associated with the stable interstellar medium 
was usually much larger than * ^ time delay x^ which is associated with the 
nebula (J. M. Rankin, private communication). This means that the effect 
of the nebular screen was small enough that the data are consistent with 
scattering in a single screen having t q ~ 10 msec at 100 MHz. 

The angular size observed at the earth due to scattering in a single 
screen with ~ 10 msec is predicted to be 0.15 arcsec at 100 MHz, ising 
Z = 2 kpc in equation V-l. This angular size scales as frequency-squared 
to 0.10 arcsec at 111 MHz. The agreement with observations is better than 
for the single-screen model only because the value of x q is smaller. The 
predicted angular size is still 50 percent larger than our observed value. 
Since our errors are only 15 percent, this is a significant difference. 

Our observations of changes in lp/lj are in qualitative agreement 
with the trends in the dispersion measure and scattering variations. 

Between the winters of 1971-2 and 1972-3 we observed that the value of 
Ip/I^, decreased from "*0. 6 to~0.4, while during this period the disper- 
sion measure showed a dramatic increase and the scattering certainly 
Increased also. Quantitatively though, the observed angular size is 
associated with a smaller interstellar time delay than the model, whereas 
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the observed value of Ip/I^ ( as in the single-screen case) predicts a 
larger time delay than the model. Since the interstellar time delay is 
not variable, this is a basic discrepancy of the observations with the 
model. Finally, the predicted decrease of fringe visibility with pulse 
phase is not observed, although the experimental uncertainty is high. 

D. Extended medium model 

From the discussion of the last two sections, it is apparent that 
either the thin screen approximation or the Gaussian irregularity spec- 
trum, or perhaps both, must be abandoned. Williamson (1972) derived 
probability density functions for the case of an extended medium in which 
the source’s radiation executes a random walk along the entire path length. 
Figure V-9 shows the impulse response as computed by Williamson. The 
rounded peak and long rise 
time are reflected in the 
predicted pulse shapes. Both 
features are due to the fact 
that with an extended medium 
there is very little "straight- 
through" radiation; the rays 
which arrive near the start of 
the pulse may come from any part 
of the scattering disk, and the 
rays which arrive from exactly 
the source direction ^ay not have 

suffered the smallest delay. This also results in an absence of fringe 
visibility dependence on pulse phase. 



Figure V-9. Impulse response of 
extended medium. 
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On this model the time delay associated with a given angular diameter 

2 

is 7r / 3 ~ 3*3 times smaller than that for the thin screen model. There- 
fore our observed angular size of 0.07 arcsec at 111 MHz corresponds 
to less than 1 msec of time delay, implying that most of the pulse broad- 
ening must occur elsewhere than in the "extended medium." Alternatively 
the observed pulse broadening of 10 to 15 msec implies an angular size 
which is much larger than that observed. 

The pulse shapes predicted using this model have rise times which 
are too long to agree with observations. The observed 111-MHz pulse from 
the Crab nebula pulsar has a steep initial rise, about 4 to 5 msec from 
baseline to peak, whereas the predicted shape using this model has a rise 
time of about 9 msec. After comparing the theoretical and observed pulse 
shapes, Williamson (1974) concludes that the extended medium approach is 
definitely excluded and that the scattering medium must be concentrated 
in a few locations along the path. 

E. Square-law structure function 

Cronyn (1974) proposed that the "structure function" of the inter- 

2 

stellar dispersion measure have the form D^(b) * D Q b where is the 
variance (mean square) of the difference in electron density between two 
lines of sight separated by a baseline b, and D q is a constant. Thus a 
Gaussian size spectrum is rejected in favor of a more realistic power 
law spectrum. It has been established that for the Interplanetary medium 
a power law size spectrum probably exists, and the same proposal for the 
interstellar medium is a natural extension of this idea. 

For a square-law structure function the group delays, due to a slower 
propagation speed through the largest irregularity sizes, are comparable 
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to the geometrical delays. The additional group delays are randomly dis- 
tributed about their mean, and they effectively randomize the dependence 
of fringe visibility on pulse epoch. 

This theory does not require the Crab nebula to be scattering screen. 
Indeed the nebula cannot be important for scattering and, in particular, 
it cannot be the cause of the variable dispersion measure and pulse shapes. 
According to Cronyn, these variations are simply manifestations of statis- 
tical fluctuations in the line-of-sight electron column density. He as- 
sumes implicitly that dispersion measure and scattering are intimately 
related as suggested by Rickett (1970). Another feature of Cronyn's 
theory is that it predicts the occurrence of dispersion measure fluctu- 
ations for all pulsars. The rms variations should depend on distance, 
relative velocities, and length of observing period. To date no changes 
other than for the Crab nebula pulsar have been observed but dispersion 
measures are being monitored at several observatories. 

Due to the group delays there is very little "straight-through" 
radiation, as in the extended medium case, and the predicted pulse shapes 
again have too-round peaks and too-slow rise times. The impulse response 
is the convolution of an exponential function (for the thin screen) with 
a Gaussian (for the random delays); it is plotted in figure V-10. The 
predicted pulse shapes cannot be 
actually compared to observations, 
however, because the essence of 
this theory is the statistical 
nature of the scattering medium. 

Therefore the predicted pulse 
shape is the average, over as 



Figure V-10. Impulse response 
for square-law structure function. 
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long a time as possible, of all the pulses which are observed. The rms 
variation In pulse shapes will continue to grow until we have observed 
the largest deviations from the average value; until a suitably large 
body of data is in hand adequate comparisons cannot be made. 

F. The Crab nebula as scattering sir 

This section addresses the question ; is the Crab nebula a scat'.ering 
screen which has a major effect on the pulsar's rudxation? That it con- 
stitutes a screen is certain because there must be some thermal electrons 
which are in some way irregularly distributed; but the screen's effective- 
ness depends on the magnitude of the density fluctuations. We will assume 
that the nebula is responsible for most of the observed pulse broadening 
and see what this implies abou the density fluctuations and scale sizes. 
Then considering what is known about the nebula we can decide whether this 
model is reasonable. 

1. Values for x^, 8^, .and scale sizes. The di? ussxon in this 

section concerns only the scattering phenomena which occ*~r in the nebula; 

the effect of the interstellar medium is considered separately. We 

adopt a value for the characteristic time delay of x^ ■ 15 msec at 

111 MHz, where the subscript N refers to the nebula. We take the path 

length through the nebula as 1 pc and let the srreen be effectively 0.5 

pc from the pulsar. Then with the total path as Z ■ 2 kpc w< have 

q ~ 4000 from equaticn V-2. Thus from equation V-l and the above values 

for t^, Z, and q, we obtain 6^ ** 2.8 milli-arcsec. The angle 6^ is 

the half-power diameter of the angular spectrum (equivalent to 0 fl of 

equation V-4). What does the value of 0., imply about structure in the 

« 

nebula? Following Cohen and Cronyn (1974) we can consider three charac- 
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teristic linear sizes in the screen which are important for scattering. 

The smallest significart size ix i. - X/20. ~ 10^ cm, the linear dis- 
tance across which one radian of phase shift occurs. Sizes smaller than 
Jt^ proc i too small a phnse shift to influence the angular broadening. 

j — 12 

The Fresnel zone size, /XZ ~ 10 cm, is the largest size which is sig- 
nificant for angular scattering; sizes large/ than fkZ cause refraction 

and have little effect on source broadening. Finally, tue scattering 

14 

disk size. Z9^ ' 10 cm, defines the distance over which we can receive 

radiation from the source. 

Such a wide range of sizes (four orders of magnitude) might be rep- 
resented by a power law spectrum of irregularity sizes in the nebula. 

The small and large lengths are equal to the gyroradii of electrons 
in a 100-pG field w*th energies of 100 and 10 Mev, re&oec lively. In 
the nebula there are electrons between these extreme energies; they would 
produce synchrotron radiation between 6 *Hz and optical frequencies. 

Could these sizes become thermal electron density irregularities? Wentzel 
(1969) suggested a way for this to happen by noting chat since streaming 
cosmic rays are scattered by their own Alfven waves, they interact racst 
strongly when they travel through one wavelength in one gyro-period. An 
investigation of the possible hydromagnetic waves in the nebular plasma 
would be interesting in that it could provide information auout possible 
density irregularity sizes. One important parameter for such a study 
would bo the thermal electron density in the nebula since this influences 
both the excitation and damping of the waves (Scheuer and Tsytovich 1970). 
The various determinations of the electron density are considered next. 

2, Thermal electron density . a. fTpansion . A firm limit can 

be set on the thermal electron density along the line of sight to the 



99 


pulsar from the lack of an observed temporal monotonic decrease in the 

dispersion measure. Since the nebula is expanding, the volume density 

systematically decreases with time. Using the known expansion rate, the 

fractional change in the dispersion measure of the nebula is calculated 

to be ~0.2% per year. Since no monotonic decrease in dispersion measure 

(i.e. column density) has been observed to a limit of about AN < 10^^ cm 

-3 

per year, the average volume density of electrons must be n^ < 1 cm 

along the line of sight to the pulsar. The contribution of the nebula to 

-3 

the total dispersion measure is then <1 pc cm if we take 1 pc as the 
path length through the nebula. This argument was first used by Drake 
(1969) to establish an upper limit to the mass of the nebula. 

b. Comparison with PSR 0525 . Other ways of determining the 
electron density involve comparisons between the Crab nebula pulsar 
(PSR 0531+21) and the nearby pulsar PSR 0525+21. They are separated on 
the sky by about 1?5 and have similar dispersion measures and rotation 
measures. The dispersion measures and rotation measures are listed in 
table V-l for both pulsars and the nebula; these values are given by 

Table V-l 


Dispersion and rotation measures 
for pulsars and nebula 

Dispersion Rotation 



PSP 0531 56.8 -42.2 

PSR 50.9 -39.6 

<1 -25.6 

Manchester (19/1). The near equality of the rotation measures of the 
two pulsars led Wilson (1974) to conclude that the radiation from 
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PSR 0531+21 is not Faraday-rotated in the nebula; rather he postulates that 
the rotation measures of both pulsars are due entirely to the interstellar 

medium. In fact he used the difference in their rotation measures (~2 rad 

-2 -3 

m ) as an upper limit to the effect of the nebula and derived n g < 0.01 cm 

along the line of sight to the pulsar within the nebula. Manchester (1971) 

notes, however, that the similarity of the two rotation measures could be 

coincidental and that a small filament in front of the pulsar in the nebula 

could account for the difference between the nebular and pulsar rotation 

measures . 

Wilson's treatment assumes that the interstellar magnetic field 

and electron density are the same along the path to both pulsars, and that 

both pulsars are at approximately the same distance. Consequently we must 

also be able to say that the difference in the dispersion measures of the 

pulsars indicates the contribution of the nebula to PSR 0531+21 *s dispersion 

-3 -3 

measure. The difference is ~6 pc cm , giving n g * 6 c* averaged through 
the nebula. But this is in conflict not only with the limit on n^ set by 
the expansion but also with the value just derived from comparing the 
rotation measures of the two pulsars. We must conclude that the inter- 
stellar medium is not the same along the two paths to the two pulsars, 
and look for other means of determining n e> 

c. Optical observations . Filaments of all sizes thread the 
nebula and it would seem that any one line of sight should encounter at 
least one. In fact there is a small filament which appears to overlie 
the pulsar. In Trimble's (1968) system it is number 159; its radial 
velocity of -162 km/sec meant, it ought to lie on our side of the pulsar 
and be about 0.2 pc in front of it. The filament is very faint, and 

_3 

probably has an electron density of about 100 cm , derived from photo- 
metric observations. It is unresolved on photographs with resolutions 
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of one arcsec and so it must be smaller than 0.01 pc. Its contribution 

-3 

to the dispersion Measure of the pulsar is thus <1 pc cm . If there 
were no other electrons along the line of sight through the nebula, the 

_3 

average electron density over the 1 pc path would be n^ < 1 cm . Com- 
bined with the limit set by the expansion, the total average electron 

•3 

density along the line of sight to the pulsar is n < 2 cm . 

€ 

For the region near the center of the nebula, Scargle (1969) 

-3 -3 

derived ~ 10 cm and suggested that the density might increase 
further from the pulsar. 

The discussion of this section has shown that the thermal electron 
density of the Crab nebula along the line of sight to the pulsar is not well 

_3 

known. However an upper limit of n ft < 2 cm can be set confidently. 

3. Filament 159 . The small filament mentioned in the previous sec- 
tion which lies between us and the pulsar could be, literally, a thin 
scattering screen. The configuration is shown in figure V-ll. In this 
section some consequences of 
the location of this filament 
are explored. We specify some 
irregularity scale sizes and 
velocities and attempt to 
justify them as physically 
reasonable. We begin by 
assuming, as stated at the 
beginning of section F, 

that most of the temporal broadening occurs in the nebula. 

a. Numerical quantities . For the fila m e n t in the location 



Figure V-ll. Filament 159 
configuration in Crab nebula. 


shown in figure V-li, we have q ~ 10000 from equation V-2. With the adopted 
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time delay of = 15 msec at 111 MHz, we calculate the apparent angular 
size of the pulsar as observed at the earth to be 0„ 35 1.1 mjLlli-.'.rc- 
sec. To proceed, we need a relation between the rms scattering angle and 
the irregularity scale sizes and density fluctuations. Derivations have 
been given by Cronyn (1970a) and Lovelace (1970) and the result is: 



d = thickness of screen. 


We have assumed that the irregularities have a typical scale size a, and 
that they are Gaussian-distributed. This may not be the most physically 
reasonable assumption, but for ease of calculations and for simplicity, 
let us continue with the Gaussian spectrum. 

We take the scale si^e as the linear distance over which one 
radian of phase shift occurs: a ~ 10^ c® (this was defined as in sec- 
tion F-l). From equation V-5 we can then calculate the rms density fluc- 

-4 -3 

tuations which could account for 0^. This yields An^ ~ 4 * 10 cm where 
we used d = 0.01 pc and the values of a and 0^ given above. If the fila- 
ment is actually thinner than 0.01 pc, the required rms density fluctu- 
ations would be larger. These calculations have ignored the regions of 
the nebula outside the filament because nothing is known of possible 
thermal electron density irregularities there. 

-3 

If the average electron density in the filament is n ~ 100 cm , 
the value of An^ Implies <0.01Z density fluctuations would be required 
to occur over a length scale of 10^ cm. Such deviations from the average 
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electron density are reasonable, and could reflect spatial variations in 
the ionization and recombination occurring in the filament. We can there- 
fore explain all of the temporal broadening as a result of density fluc- 
tuations in the filament. 

b. Time delay varia ti ons . Another requirement of the model is 

that changes in x^ which are impl: ed by our observations of changes in 

the ratio Ip/I^ should also follow from physically reasonable variations. 

From the graph of figure V-4, the change in the value of Ip/l^ from 0.6 

(observed in 1971-2) to 0.4 (observed in 1972-3) implies a change in x^ 

from 15 to 25 msec over a time interval of one year (3 x 10^ sec). This 

in turn yields a change in 0^ of a factor of ~1.3, and a corresponding 

—4 —3 

variation of An from 4 to 6 * 10 cm . 
e 

The pulsar's proper motion is ~100 km/s (Trimble 1968) and the 
filament also has a proper motion of ~30 km/s (Trimble 1968) in a direction 
approximately perpendicular to that of the pulsar. It is not unreasonable 
to suppose that relative motion of the pulsar, filament, and observer 
could carry the line of sight through regions of the filament where the 
ionization is more or less patchy. These quite minor variations in the 
mean electron density of the filament can easily account for the observed 
variations in scattering parameters. 

c. Changes in dispersion measure . Can irregularities in the 
filament also be used to explain the observed variations in dispersion? 

The variations reported by Rankin and Counselman (1973) can be viewed as 
random, overlapping "events" which result in increases and decreases in 
the column density of thermal electrons. These changes could conceivably 
be due to "clouds" of appropriate size, velocity, and thermal electron 
density moving across the line of sight. Rough values of the required 
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cloud parameters were estimated by Rankin and Roberts (1971) based on 
"typical" magnitudes and durations for the "events." Subsequent obser- 
vations revealed that the variations are much more complicated, as shown 
in the two-year span of data presented by Rankin and Counselman (1973). 

It is possible that within the filament there exist large-scale 
irregularities which are responsible for the changes in dispersion. These 
irregularities would be of the order of ,or larger than ,the linear diameter 
of the scattering disk in size. Smaller sizes will be numerous enough that 
as many will move into as out of the area of the disk and no changes in 
dispersion measure due to these would be observed. Let us suppose there 

are clouds of size a^ — 2 a.u. in the filament, and let each have an aver- 

-3 

age density of n^ ^ 200 cm . Then with a 50Z filling factor we can pro- 
duce the observed emission measure of 200 pc cm ^ and the limit on the 

-3 

average density of n^ ~ 100 cm . If the clouds have velocities of 

v ~ 100 km/s (which could be due to the proper motion of the pulsar and 

filament), one "cloud" would paes across the scattering disk in a time 

t - a^/v * 3 x 10^ sec ~1 month, and its passage would cause a change in 

15 -2 

electron column density of AN * n^a^ “ $ x 10 cm . The observed dis- 
persion variations have approximately these time scales and magnitudes. 

Until more extensive observations are available so that statis- 
tical methods can be used, these simple calculations must remain suspect 
because there are too many unknown and arbitrary parameters. 

d. Scintillations . Temporal variations in the pulsar's in- 
tensity are also predicted due to irregularities in the filament. A dif- 
fraction pattern would be formed by interference between various 
scattered rays. The scale of the pattern will be on the order of 
l . ~ 10*® cm. The scintillations would esdtibit a time scale of 
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3 

i,/v ~ 10 sec and be correlated over a bandwidth of 1/2itt„ ~ 10 Hz at 
<p N 

111 MHz. The observable variations will be discussed In the next sec- 
tion. 

e. Interstellar scattering . On the model described above, the 
source which appears at the edge of the nebula at 111 MHz has been broaden- 
ed Into an apparent angular size which, viewed from the earth, is a 
few milli-seconds of arc; the pulse shape has been broadened to display 
a characteristic time delay of *-15 msec; and intensity scintillations have 
been built up in the radiation. 

The source is now larger than the critical size required for 
interstellar intensity scintillations (Cohen and Cronyn 1974), though its 
passage through the interstellar medium will not alter the already-present 
variations. Interstellar scattering will additionally broaden the angular 
size and pulse shape of the source, however. If we adopt Williamson's 
model for an extended medium (described in section D) we can achieve an 
angular diameter of 0.07 arcsec while convolving only 0.7 msec of temporal 
broadening into the pulse shape. An earth-based observer would then per- 
ceive a source which has the requisite combination of temporal and angular 
broadening. 

On this model, or any model in which appreciable temporal broad- 
ening occurs in the nebula, narrow-band structure observed in the pulsar's 
spectrum will be indicative only of the nebular time delay. This model 
therefore predicts that the total observed time delay (t q ) and the width 
of the frequency structure (fjj) will not be in the relation 2irT o f^ • 1 
as was observed for the Vela pulsar by Backer (1974). 

At 111 MHz, where ~ 15 msec, the decorrelation bandwidth is 
only f^ ^ 10 Hz. At higher frequencies f^ approaches a feasible observing 
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bandwidth and this prediction is testable* Hugu nin et al. (1969) observed 
this pulsar at 151 MHz and observed no frequency structure in the radiation. 
However they used a passband of 0.1 MHz which would have blurred any scin- 
tillations. 

G. Summary of models 

The five models presented above all attempt to explain the observed 
features in the radio radiation from the Crab nebula pulsar. The major 
observation which any model must explain is the combination of angular 
and temporal broadening which have both been observed at 111 MHz. Based 
on the temporal broadening, the simplest model, a single thin screen, 
predicts an angular size which is about a factor of two larger than the 
size measured with VLBI. A similar discrepancy exists for the two-screen 
model. Also for both models, the predicted decrease of fringe visibility 
with pulse phase is not observed. The extended medium model of Williamson 
again predicts an angular size which is much too large to agree with the 
VLBI observations. To evaluate the structure function of Cronyn, more 
observations of other pulsars are still needed. Finally, the model in 
which nearly all of the temporal broadening occurs in the nebula and only 
a few milliseconds are added in the interstellar medium does predict the 
observed combination of temporal and angular broadening. 



CHAPTER VI 


OBSERVATIONS OF OTHER PULSARS 

This chapter discusses the observations of pulsars other than the 
Crab nebula pulsar* The first part of the chapter is a presentation of 
the results and their implications* The second part is a discussion of 
various interstellar scattering effects. 

A. Results and implications 

1. Strong pulsars . a. Observations and processing . Four strong, 
low-dispersion pulsars were observed during each session; table VI-1 lists 
the pulsars' parameters. Although many tapes on each pulsar were recorded 


Table VI-1 

Low-dispersion pulsar parameters 


PSR 

Period 

(sec) 

Equivalent 
width (msec) 

Duty cycle 
% 

Dispersion 
measure (pc cm 

Gate 

width (msec! 

0834+06 

1.274 

20 

4 

12.9 

50 

0950+-08 

0.253 

10 

4 

3.0 

10 

1133+16 

1.188 

20 

4 

4.8 

50 

1919+21 

1.337 

25 

2 

12.4 

25 


during every observing session, results are reported here for only a few 
runs. One reason for the small number of runs is that since all pulsars 
are highly iriable sources, many times they were too weak to be detected 
in a three-minute integration. If the correlated signal could be detected, 
individual pulses were usually quite easily visible on the total power 
records. Due to the complexity of having three stations on the air 
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simultaneously, it was not possible to revise the schedule when it appeared 
that a pulsar happened to be weak at the time we observed. 

The other reason for the few runs is that we consider here only 
those observations for which fluxes on all three baselines are available. 
This is again because pulsars are highly variable sources, and the fringe 
visibility on the long baseline can be calculated only if a simultaneous 
measurement of the total flux was made using the short baseline. This 
excludes the 1972 January, March, and December data when only one base- 
line was active. As in the case of the Crab nebula pulsar, the fringe 
visibility, y, is defined as the ratio of long-baseline (S p ^) to short- 
baseline (Sp) flux, y is also the fringe visibility of the pulse energy, 
since the gate width was the same for all baselines. Finally, only the 
111-MHz observations of pulsars are reported here. At 196 MHz few three- 
baseline observations are available, and at 74 MHz the spectra of most of 
the pulsars have turned over making the pulsars too weak to detect. 

The gating method was used for all the pulsars listed in table 
VI-1. If, through an error or mishap, the pulse epoch was not determined, 
searches for the correct epoch were not made because of the large amount 
of computer time involved. The gating method was necessary on every pulsar 
observation because the enhancement of the signal-to-noise ratio meant the 
difference between detection and non-detection. Several runs were proces- 
sed both without a gate (treating the pulsar as a continuum source) and with 
a gate, and only those runs processed with the gate showed fringes. Table 
VI-2 gives the detection limits at 111 MHz for each pulsar, both for a 
full-tape integration and for a single pulse. 

b. Average fluxes and fringe visibilities . Table VI-3 lists 
the average fluxes, fringe visibilities, and angular sizes at 111 MHz for 
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Table VI-2 

System sensitivity (5o level) for pulsar 
at 111 MHz, by baseline 



0834+06 

141 

380 

1.4 

2.6 

3.7 

17 

31 

44 

0950+08 

711 

470 

1.6 

2.9 

4.2 

42 

76 

112 

1133+16 

151 

540 

1.7 

3.0 

4.4 

20 

37 

55 

1919+21 

134 

1500 

8.9 

8.9 

14.6 

61 

103 

170 




Average pulsar parameters. 111 Milz 
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strongest pulses could be use-, there were typically a dozen or so pulses 
in each run. 

It was found that there i3 no difference between the visibilities 
deri\ *d from the single pulses and those derived from the average fluxes. 
This is in contrast with the predictions of Cronyn (1972) that thf fringe 
visibility should decrease for longer coherent integration times, and in- 
dicates that decorrelation across the long baseline due to interplanetary 
scintillations was not as severe as indicated by the rms phase deviations 
presented in table II-7. None of the pulsars lies in the ecliptic plane, 
the closest being PSR 0950+08 which is 5° above the plane, so that a de- 
crease of scintillations away from the plane might explain the lack of an 
appreciab > effect on these pulsars. Also, most of the observations were 
made at solar elongations greater then 90° which could account for a reduced 
solar wind effect. 

d. Single-pulse amplitude and phase plots . All of the pulsars 
we observed have pulse widths less than 0.2 sec (the length of one record) 
so that the correlation amplitude of a single record reflects the pulse 
intensity. To show that the same amplitudes were observed on each base- 
line, the single-pulse fringe amplitudes were plotted as a function of 
time or pulse number. The amplitudes observed on all three baselines 
were superimposed and the pulse-energy measurements made at Arecibo were 
also plotted when they were available. Figure VI-1 shows the pulse-to- 
pulse amplitudes and phases for an observation of PSR 1919+21 in 1971 
December. The data were smoothed by a five-pulse running average. The 
dramatic increase in flux during the latter half of the run illustrates 
the variability of the pulsar flux and also shows that both the unresolved 



LONG BASELINES (ARECIBO-NRAO aARECIBO SUGAR GROVE ) 
SHORT BASELINES (NRAO- SUGAR GROVE) 

ARECIBO 



Figure VI-1. Pulse-to-pulse amplitude ana phase at 111 MHz for ’SR 1919+21 iu 
1971 December. The data were smoothed with a five— pulse running average. The 
total length of the pulse series is three minutes. 
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and partially-resolved fluxes exhibit the same variations. The pulse am- 
plitudes plotted in figure VI-2 offer a bet «sr example of the closeness 
of the variations on all baselines. This data was obtained for PSR 0834406 
in 1973 February; although Arecibo pulse-energy measurem its were recorded 
they are not yet available for comparison with the interferometer data. 

The amplitudes and phases were smoothed with a five-pulse-wide filter to 
narrow the noise spectrum. 

The presence of very strong pulses distorts the distribution of 
fringe amplitudes: pulsars display amplitude distributions which have 

long high-amplitude tails in contrast with the distributions for continuum 
sources which are Kicean. Figure VT-3 is a plot of the distributions of 
fringe amplitude and phase for the quasar 3C286 and for the pulsar PSR 0950408. 
The amplitude distribution for PSR 0950408 is similar to other histograms pub- 
lished for this pulsar (Smith 1973). The two phase distributions are both 
approximately Gaussian. This illustrates that the correlated pulse ampli- 
tudes have the same characteristics as those observed with a single tele- 
scope. 

e. Polarization . Most pulsars are highly polarized, so that 
even though our instrumental polarization was very high, we could measure 
the linear polarization of pulsars. During the 1973 February experiment 
we made several observations in linear polarization on PSR 0834406 and 
0950408. The ^ luxes were corrected for instrumental polarization accord- 
ing to the equations derived in Appendix D. We observed less than 15Z 
linear polarization for PSR 0834406, and 90±20Z linear polarization for 
PSR 0950408; these numbers are derived from the average of the flux ob- 
served during three runs on each pulsar. 



SHORT BASELINE ( NRAO- SUGAR GROVE ) 

PSR 0834 +06 LONG BASELINES (ARECIBO-NRAO & ARECIBO* 

25 FEBRUARY I9V3 SUGAR GROVE ) 
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Figure VI-2. Pulse-to-pulse amplitude and phase at 111 MHz for PSR 0834+06 In 
1973 February. The data were smoothed with a five-pulse running average. The 
total length of the pulse series Is three minutes. 
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Figure VI-3. Distribution of amplitude and phase for PSR 0950H)8 
and 3C286. Both runs were made at 111 MHz in 1973 February. There 
are a total of 710 pulses for the pulsar and 900 samples for the 
quasar. The amplitude scale is in units of the mean amplitude. 

The phase distribution is wider for the pulsar because the signal- 
to-noise ratio is lower (see equation 11-22, page 41). 
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Manchester et al . (1973) collected polarization measurements for 
20 pulsars at many frequencies and concluded that most pulsars exhibit a 
constant fractional linear polarization below some critical frequency 
while the polarization decrease above that frequency. Although our measure- 
ment for PSR 0950+08 is consistent with the limiting polarization of 65Z 
assigned by Manchester et al * , it may be that this pulsar does not reach 
the limiting polarization until below 100 MHz. Our measurement for 
PSR 0834406 is not consistent with the measurements of Manchester et al . 

They find ~25Z polarization at 151 MHz and believe that the limiting 
polarization has not been reached. According to this, we should have 
observed about 25Z polarization; however, the high instrumental polar* 
ization precluded more accurate measurements* 

2* Weak pulsars . We hoped to observe substantial resolution and 
to measure an apparent angular size for distant pulsars other than the 
Crab nebula pulsar to establish a dependence of angular size on distance. 

For this reason, ve observed four weak, high-dispersion pulsars at 111 
and 74 MHz during 1972 December and 1973 February. Table VI-4 lists the 
pulsars 1 parameters, as given by Lang (1971c). These pulsars were pro- 
cessed as continuum sources because they all have little flux which is 
pulsing at low frequencies. No correlated flux was detected for any of 
the pulsars on either the long or short baseline. 

Table VI-4 

High-dispersion pulsar parameters 


PSR 

Period 

(sec) 

Dispersion measure 
(pc cm"3) 

l 

b 

1858+03 

0.655 

402 

3712 

lo!? 

1933+16 

0.359 

158 

52.4 

-2.1 

1946+35 

0.717 

129 

70.6 

5.0 

2003+31 

2.111 

225 

69.0 
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There are many reasons for not observing fringes, but ruling out 
errors in recording or processing leaves only two possible reasons. First, 
the source size may be comparable to the fringe spacing and therefore the 
source is partially resolved by the interferometer. This possibility could 
obtain for these four pulsars since their asymmetrically broadened pulse 
shapes (Lang 1971c) indicate that interstellar scattering should also pro- 
duce a broad angular size. Second, the source may be actually smaller than 
the fringe size but too weak to detect. This could occur if the emitted 
spectra of these pulsars turn over at frequencies above about 100 MHz; then 
the flux probably would have been below our detection limits on all base- 
lines. 

Seiber (1973) compiled all the available flux measurements for PSR 
1933+16 and published a spectrum which exhibits a turnover near 300 MHz. 
Therefore the flux was well below our detection limit on all baselines. 
Insufficient flux measurements have been published for the other three 
pulsars to enable this same statement to be made, but some limits can be 
established. Recently McLean (1973) published spectral indices for PSR 
1858403 and 1946+35 using observations at 922 and 408 MHz. Extrapolating 
to our frequencies we obtain the fluxes listed in table Vl-5. Our upper 
limits on the flux indicate that the spectra of these two pulsars could 
turn over below 408 MHz. Alternatively, the pulsars could have been weak 
at the times we observed. 

Table VI-5 

Fluxes of PSR 1858+03 and 1946+35 


PSR 

Spectral 

Index 

Frequency 

(MHz) 

Extrapolated 
Flux (Jy) _ 

Observed 
Flux (Jy) 

No. of 
runs 

1858403 

-2.9 

111 

3.3 

<4±2 

4 



74 

10.6 

<10±5 

2 

1946+35 

-2.7 

111 

4.0 

<1.9±1 

8 



74 

11.9 

<4.5±2 

5 
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Without better knowledge of the spectra we cannot interpret the lack 
of detection on the long baselines quantitatively. The apparent angular 
sizes of these pulsars are probably partially resolved but actual total 
flux measurements rather than limits are needed to calculate even a limit 
on the fringe visibility. 

B. Discussion of interstellar scattering effects 

The plane wave which starts out from the pulsar and propagates toward 
the observer is perturbed during itc passage through the irregular inter- 
stellar medium. Both amplitude and phase fluctuations are imposed on the 
wave, and each manifests itself in a different way. Amplitude variations 
(due to interference effects) are the basis for the observed scintillation 
phenomena, and phase perturbations in the wave give rise to predicted, but 
not yet observed, apparent source mot: on. The multipath scattering also 
produces angular and temporal broadening of the radiation; this aspect was 
discussed in Chapter V. Each of these phenomena will be discussed with 
regard to the pulsar observations. 

1. Scintillations . The characteristic time delay discussed in 
Chapter V also implies a characteristic frequency width, » f^, 

over which intensity fluctuations are correlated. The quantity f^ is 
called the "decorrelation bandwidth," and in order for intensity scintil- 
lations to be observable the instrumental bandwidth must be smaller than 
fp so that the frequency-dependent diffraction pattern is not blurred or 
"decorr elated." At 111 MHz, estimates for f^ are 1 to 2 kHz for all the 
low-dispersion pulsars we observed except PSR 0950+08 which has f ~ 100 kHz. 
The Mark I bandwidth of 350 kHz is thus many times wider than f^ and for 


this reason it is unlikely that the observed variability is due to 
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interstellar scintillations. Interplanetary scintillations are also 
probably not the cause of the flux variations either, because the tine 
scales for IPS are very short. Instead, we attribute all flux changes 
to intrinsic source variability. Possibly for PSR 0950408 some variations 
are due to scintillations since the observing and scintillation bandwidths 
were quite close. 

The observed fluctuations in intensity occurred over time scales from 
s a pulses to over one year. The similarity between the individual 
pulse energies measured at Arecibo and the single-pulse correlated fluxes 
(figure Vl-i) implies that the variability observed with the interferom- 
eter reflects the true pulse-to-pulse fluctuation. Over monthly intervals 
the flux changes are also probably intrinsic to the source. Huguenin et al . 
(1973) note that changes of up to a factor of 10 occur on characteristic 
time scales which are different for each pulsar. Over several months we 
observed large changes in flux, but from day-to-day within an observing 
session the average flux stayed approximately the same. 

2. Source motion . The possibility of apparent position variations 
due to interstellar scattering is discussed by Lovelace (1970). This 
phenomenon is analagous to the effect called "seeing" by optical obser- 
vers. Lovelace notes that the apparent position changes should usually 
be caused by the largest irregularity sizes in the medium. For a Gaussian 
irregularity spectrum there is no predicted apparent source motion because 
only for a steep power law spectrum (with spectral index >2) are there 
sufficiently large irregularity sizes to produce the effect. It is dif- 
ficult to predict the magnitude and time scale of the variations without 
knowledge of the Irregularity size spectrum. We would expect, however. 
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that observable position variations would be at least several times larger 

than the scattering size, 6 , and that the time scale would depend on the 

s 

linear size and relative transverse velocity of the irregularity which 
produces the effects. Typically, we expect an excursion of one scattering 
size in one intensity scintillation time. 

Since our interferometer fringes are on the order of the apparent 
angular size of the pulsars, we would observe phase excursions of at least 
several turns over time intervals of up to several years depending on the 
physical parameters of the irregularities. Without any more precise 
specifications, this effect of interstellar scattering is hopelessly en- 
tangled with ionospheric refraction effects. Over our long baseline (2500 km) 
ionospheric scintillations are completely decorrelated and at 111 MHz we 
could observe changes in the source's apparent direction of several hun- 
dred fringes (equivalent to about one minute of arc) due to the ionosphere. 
Such deviations are within the interferometer's response pattern but changes 
this large could completely mask source position fluctuations which are due 
to interstellar scattering. 

Our observations of pulsars show no position variations which cannot 
be attributed to the ionosphere. To properly investigate source motion, 
differential techniques should be used in which one measures the phase 
difference between the pulsar and another source also in the telescope 
beam. Such a technique was successfully used to put a limit on the posi- 
tion difference between the Crab nebula pulsar and compact source (Chap- 
ter IV, section C-4). 

3. Angular broadening . Apparent angular sizes at 111 MHz for a 
total of five pulsars were measured during the experiment. If we look at 
equation V-5 (page 103 » we see that (apart from wavlength) the angular 
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size depends directly on the source distance (Z) and the ros electron 
density variations (An e ) , and inversely on the irregularity size (a) . If 
we assume that the irregularities are associated with the electrons res- 
ponsible for the dispersion (D) and that the interstellar medium is similar 
along lines of sight to all pulsars, then the observed angular sizes should 
scale as the square root of the dispersion. The observed angular sizes at 
111 MHz are plotted against dispersion measure in figure VI-4, and the 

values for 9 , D, and Z are listed in table VI-6. The distances are those 
s 

assigned by Prentice and ter Haar (1969) who took account of the known HII 

Table VI-6 

Pulsar sizes (at 111 MHz) and distances 


PSR 

9 

(arcsec) 

_3 

D(pc cm ) 

Z (pc) 

0531+21 

0.07 ±0.01 

56 

2000 

0834+06 

0.06±0.01 

13 

400 

0950+08 

0.05±0.03 

3 

60 

1133+16 

0.06±0.04 

5 

130 

1913+21 

0.06±0.03 

13 

250 


regions and early-type stars along the line of sight to the pulsars. In 
figure Vl-4 a line of slope one-half is drawn for both 111 and 74 MHz. 
The line for 111 MHz passes through the point for the Crab nebula pulsar 
since this is the best-determined size. The nearby pulsars appear to be 
too large at 111 MHz, but the large error bars do Include the predicted 
values. 

What do the changing angular sizes (table VI-3) indicate? A de- 
crease in angular size means that the rms density fluctuations along the 
line of sight have smoothed out, or that larger scale sizes drifted into 
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at 111 MHz. Lines have slope one-half, and at 111 MHz the line passes through 
the size of the Crab nebula pulsar. The 74-MHz line is shown in order to in- 
dicate the wavelength dependence. The 1973 February points from table VI-3 
are plotted. 
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the line of sight. For the pulsar observations in 1973 February, D. C. 
Backer and G. A. Zeissig recorded individual pulse intensities simultan- 
eously with our experiment at Arecibo and NRAO so that the drifting dif- 
fraction patterns could be observed. The data have not yet been analy .ed, 
but it will be interesting to see if the patterns are different during the 
two days when different values for the angular size were measured. We 
predict that the observed patterns should be larger on the day when the 
pulsar was smaller. 

The question remains, what in the interstellar medium causes the 

scattering? The evidence from the angular broadening is that there are 

q 

myriad irregularities with linear sizes ~ 10 cm, the distance over 
which one radian of phase shift occurs. The latest diffraction pattern 
measurements (Slee et al. 1974) indicate a similar irregularity size. 

Such sizes are much too small to be observed with the conventional tech- 
niques used for studies of the interstellar medium. The existence and 
prevalence of small scale sizes are strongly indicated by various obser- 
vations of pulsars using quite disparate techniques, and new theories 
concerning the interstellar plasma need to consider this fact. 



CHAPTER VII 


SUMMARY, CONCLUSIONS, AND SUGGESTIONS FOR FURTHER WORK 

This thesis presented angular size measurements for five pulsars, 

Wj.th the body of data heavily weighted toward the Crab nebula pulsar. A 
visibility function for< the Crab nebula pulsar which is consistent with a 
Gaussian brightness distribution was measured in an early 144-MHz experi- 
ment. Assuming a Gaussian distribution at all frequencies, we have also 
made estimates of apparent angular sizes at 196, 111, and 74 MHz. These 
sizes scale as wavelength-squared in accordance with interstellar scatter- 
ing theories, and the angular sizes for all five pulsars show a dependence 
on the square root of distance (or dispersion measure) which is consistent 
with theoretical predictions. 

We have also measured the flux of the compact source in the Crab neb- 
ula at 196, 111, and 74 MHz, and find that the fluxes lie on the extrapo- 
lation from high frequencies of the pulsing flux. Using differential fringe 
phase measurements we have tightened the limit on the coincidence of the 
pulsar and compact source to a separation of ~2 a.u. at the distance of 
the nebula. 

We have shown that the apparent angular size and temporal broadening 
of the Crab nebula pulsar cannot both be explained with a single thin 
screen model for interstellar scattering. A two-screen model improves 
the situation, but equal-strength screens predict results inconsistent 
with observations. It appears that for this pulsar, the observed com- 
bination of angular and temporal broadening may be intimately related to 
the presence of the nebula. A model in which the small filament lying in 
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front of the pulsar Is responsible for the temporal u^or.-j^ning while the 
extended Interstellar medium causes the angular V:o \ler„i <p. does enable 
the otherwise-discordant angular and temporal o be explained. 

Further observations of this pulsar, and comparisons t.o observations of 
other pulsars, would help to determine if the Crab <?.-'bula does indeed have 
a unique effect. 

The results presented in the previous chapters naturally lead us to 
ask many more questions about the Crcb nebula, about interstellar scatter- 
ing effects, and about interplanetary scintillations. Possible future 
projects involving the study of these media include the following. 

1. The theories which propose the Crab nebula as a scattering screen 
predict that the narrow-band frequency structure in the radiation of the 
pulsar will not follow the relation 2trxf ** 1. This prediction could be 
tested with observations made near 400 MHz with a frequency resolution of 

3 kHz or less. If scintillations are observed which are not in accord 
with the time delay then we could conclude that the Crab nebula does play 
an important role as a scattering screen. 

2. The total flux of the pulsar (i.e. the compact source) at 111 MHz 
does not appear to vary as substantially as the pulsing flux, although we 
were able to observe only once per month. A program to monitor the com- 
pact source over short time intervals with an interferometer could de, 
mine if the variations in total flux and pulsing flux are correlated. This 
would help determine how much of the variability is attributable to scat- 
tering and how much Is intrinsic. 

3. The fhsence of small-scale structure in the Crab nebula other 
than the pulsar was inferred through arguments concerning our flat visi- 
bility curves, the total flux agreement between pulsar and compact source 
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at 196 MHz, and the position agreement between pulsar and compact source 
at .*11 MHz, More complete coverage of the (u,v) plane at these frequencies 
would pro\id; a good test of our inferences, and also of our assumption 
that the brightness distribution of the scattering disk is Gauss5an. A 
VLBI experiment employing fully-steerable telescopes and baselines at 
least as long as those used in this project would he needed. 

Another way of investigating the small-scale structure in the nebula 
would be to determine the spectra of the radiation from -he filaments by 
obtaining high-resolution synthesis maps of the nebula at several fre- 
quencies between 1 and 10 GHz. Such a project would be an extension of 
Wilson's (1972) work to higher resolutions and more frequencies and could 
be done using the 35-km interferometer at NRAO or the proposed Very Large 
Array (VIA) . 

4. The dispersion measures of many pulsars are being monitored at 
several observatories. Detect : a of variations in this quantity would 
support Cronyn’s (1973) theory of a power law spectrum of irregularity 
sizes. If variations are detected, moi itoring of otl ar pulsars' disper- 
sion measures could provide more statistical information on the dynamics 
of large irregularity sizes in ^he interstellar medium. Also, a detection 
of dispersion measure variations in pulsars other than V Crab nebula 
pulsar would alsd indicate a neea for re-examination of the models which 
propose that the nebula is responsible for the variations. 

5. A concentrated program to perform simultaneous observations of 
drifting diffraction patterns and angular sizes for many pulsars would 
yield much-needed statistical daia on th* smaller irregularity sizes. 

Th ^ observations presented here have shown that interstellar angular 
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scattering is an observable, and possibly changing, phenomenon. There is 
a need for simultaneous observations of every possible observable parameter 
of pulsars because of their great variability. 

6. An investigation of the angular extent, time scale, and wavelength 
dependence of predicted apparent source motion could be done using an inter- 
ferometer and differential phase techniques. This also would provide in- 
formation on the sizes and velocities of the largest irregularities in the 
interstellar medium. 

One should ask hew or if the experiments and techniques used for the 
present work could be improved upon to facilitate future investigations. 
Unfortunately, a consideration of limiting instrumental effects shows that 
no order-of magnitude increase in sensitivity or accuracy can be made in 
the study of pulsars 1 angular structure for the following reasons. 

To measure the apparent angular diameters more accurately, we can in- 
crease the interferometer’s sensitivity or its angular resolution. For 
better sensitivity we require a wider bandwidth, a longer integration time, 
or a reduced system temperature. None of these can be achieved at low fre- 
quencies, however. Wide bandwidths cannot be used because of interference 
problems; the integration time cannot be increased because of ionospheric 
and interplanetary phase fluctuations which decorrelate the signal; and 
the syFtem temperature carnot be reduced because it is dominated by the 
galactic background at frequencies below about 200 MHz. Even worse, for 
observations of the Crab nebula pulsar the nebula dominates the system 
temperature at all frequencies. Conventional means of increasing sensi- 
tivity do not work at low frequencies for interferometry. 
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Instead of reducing the noise level, we could perhaps increase the 

resolution so that the relative error in fringe visibility is smaller. 

In raising the observing frequency, ve must realize that the interstellar 

2 

angular scattering effect decreases as X while the fringe spacing of an 
interferometer only increases as X/D where D is the baseline length. Thus, 
the relevant parameter for interferometric observations of interstellar 
scattering effects is XD. (The dependence of angular scattering on dis- 
tance is also important but the effect is smaller.) To illustrate the 
numbers involved, table VII-1 lists the parameters for our observations 

Table VII-1 

Observing parameters for Crab nebula pulsar 


2 

Freq. Wavelength Y XD(m ) 


196 MHz 

1.5 m 

0.85 

3.8 

x io 6 

111 MHz 

2.7 a 

0.68 

6.9 

x 10 6 

74 MHz 

4.0 a 

0.27 

10.2 

x io 6 

of the Crab nebula pulsar. 

From this 

table we see that 

to avoid complete 

resolution of a pulsar, we 

require XD 

7 2 

& 10 m and to achieve partial 


6 2 

resolution we need XD £ 10 m . At a frequency near 400 MHz, a XD product 

6 2 

of 5 x 10 m requires a baseline D ~ 6700 km. Any baseline is limited to 
the diameter of the earth, and more reasonably to an earth radius so that 
mutual visibility of sources is possible. A fundamental limit to the base- 
line is therefore D & 6400 km. We can hope to observe substantial reso- 
lution of the apparent sizes only of the most distant pulsars. 

The possible projects mentioned above which require VLBI techniques 
could be done using existing telescopes which have low-frequency capabilities. 
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Table VI I- 2 lists a few possible baselines using such telescopes at 200 MHz. 
The advantage of having fully-steerable dishes to provide good (u,v)-plane 


Table VII-2 

Possible baselines at 200 MHz 


Baseline length 



coverage indicates that the Jodrell Bark to NRAO baseline nay be the best 
for the experiments proposed. When more definite programs are designed, 
the sensitivity of each baseline should also be taken into account. 

The work presented in this thesis has only begun to untangle the many 
effects imposed on pulsars' radiation by the irregular interstellar medium. 
There are many more inquiries to be formulated and investigations tc be 
undertaken which will shed light ou the structure and dynamics of the In- 


terstellar medium 



APPENDIX A 


SCHEDULING AND PREPARATION PROCEDURES 

A, Pre-scheduling program 

Obtaining observing time at more than one station simultaneously re- 
quires cooperation and coordination between the experimenters and telescope 
schedulers. Once time has been allocated, the experimenters must schedule 
observations such that all stations can point their telescopes to the same 
sources at the same time, and that there is ample time allowed to move each 
telescope from source to source. 

We use a computer program called PRESKED which has as input a list 
of sources, station locations, and telescope limit stops. It calculates 
the times at which each source rises and sets at the selected stations. 
"Rise* 1 and "set” define the time interval when the source is within the 
telescope’s limits. The mutual visibility of a source at all the stations 
is also calculated. The program output is a list of the rising and setting 
times for each source at each telescope, accompanied by a bar graph where 
the mutual visibility time is blocked out in UT. This facilitates sched- 
uling individual observations by showing at a glance which sources are "up" 
at which times. 

B. Scheduling program 

We employ computer programs for scheduling purposes also. The main 
program e xami nes a proposed observing schedule and checks each run against 
each telescope’s hour angle and elevation limits and sl^v es. It flags 

runs which cannot be made and also indicates the error, election 
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limit or azimuth slew rate. Based on this program* s output, a new schedule 
is tried. After several iterations a schedule acceptable at all stations 
is finalized. The program produces a copy of an observing schedule for each 
station, which lists for each observation the starting time and source posi- 
tion and other information. The same program prints a label for each tape 
which is scheduled to be recorded and also punches the input cards for the 
processing program. Because communication between the stations is not always 
possible, changes to the schedule are rarely made during the course of the 
experiment. Notes are added to the schedule during the experiment, such 
as which tapes started late, interference levels, and clock drift rates. 

The annotated schedule then serves as a station log. 

C. Station set-up 

Each station in the experiment has its own peculiar set-up which de- 
pends on available equipment and normal operating procedures at that tele- 
scope. Figure A-l is a block diagram of a typical station which i3 set up 
for pulsar VLBI Mark I observations. Usually a 30-MHz IF signal is avail- 
able from the front end of the telescope and is delivered to the Mark I 
controller. Extra equipment is necessary for pulsar and total power record- 
ings. Pulse epochs are measured (several millisecond accuracy) at one or 
more stations and individual pulse intensities are recorded simultaneously 
with the Mark I recording. The total power in the Mark I bandwidth is also 
monitored and system calibration observations are made periodically to 
measure receiver temperature and dish efficiency. Table A-l lists the 
adopted system parameters for each observing session. 
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Table A-l 

Adopted system parameters for low-frequency experiments 


Station 

Parameter 

Freq 

Nov 71 

Dec 71 

Jan 72 

Feb 72 

Mar 72 

Dec 72 

Feb 73 


Receiver 

196 

i 500 

280 

300 

200 

350 


. 


temp . (K) 

111 

700 

650 

600 

600 

500 

400 

400 



74 

— 

— 

— 

— 

— 

600 

600 

Arecibo 

Efficiency 

196 

1.3 

1.3 

1.3 

1.3 

1.3 

- 



(K/Jy) 

111 

3.5 

3.5 

3.5 

3.5 

3.5 

2.6 

2.6 


1 

74 

— 

— 

— 

— 

— 

4.6 

4.6 


Receiver 

196 

350 

350 

350 

350 

350 





temp. (K) 

111 

170 

170 

170 

170 

170 

200 

200 



74 

— 

— 

— 

— 

— 

400 

400 

NRAO 

Efficiency 

196 

.93 

.93 

.93 

.93 

.93 



T 


(K/Jy) 

111 

.96 

.96 

.96 

.96 

.96 

1.0 

1.0 



74 

— 

— 

— 

— 

— 

.60 

.60 


Receiver 

,Q 6 

375 

375 

375 

375 

375 

— 

— 


temp. (K) 

lxi 

675 

675 

675 

675 

675 

930 

530 

Sugar 


74 

— 

— 

— 

— 

— 

900 

610 

Grove 

Efficiency 

196 

.31 

.31 

.31 

.31 

.31 


— 


(K/Jy) 

111 

.30 

.30 

.30 

.30 

.30 

.35 

.35 



74 

— 

— 

— 

— 

- — 

.35 

.37 
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The receiver configuration and observing techniques varied substan- 
tially from one experiment to the next. As we gradually learned better 
ways of handling the many potential problems involved in low-frequency 
VLBI observations, we altered the original receivers and other equipment 
and made changes in calibration procedures. Some of the items shown in 
figure A-l should be mentioned as examples of the special circumstances 
encountered. 

The feeds at NRAO and Arecibo were the standard low-frequency feeds; 
at Sugar Grove a special 74-MHz feed was built by A. Youmans for our ex- 
periments. At low frequencies the size of the feed can be quite large and 
interaction between parts of the feed and feed-support legs can result in 
an unusual antenna response. This is probably the cause of the asymmetric 
74-MHz efficiency curve at NRAO (see figure II-8, page 33). 

Tunable RF filters were used to reject as many interfering signals 
as possible before mixing. An interference survey was done at each station 
in the fall of 1972 to identify the exact band near our nominal frequencies 
which had the least interference. As a consequence of interference signals 
which did get through, the clipper would often saturate making some data 
meaningless. Big, 250-mLlli-watt IF amplifiers were obtained for use at 
each station to insure linearity of the system even when strong interference 
was present. 

The galactic background temperature is very irregular over the sky, and 
since it forms a substantial fraction of the total system temperature at low 
frequencies, measurement of the total system emperature during each run is 
very important. Calibration of the last two experiments was done using a 
variation of the "Y-f actor” method. To determine the system temperature. 
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we measured the amount of attenuation required tc keep the signal level 
constant. For the earlier experiments, system temperatures were estimated 
using information gained from the later experiments. 

At 100 MHz, the natural fringe rate on the short baseline (NRAO to 
Sugar Grove) is only about 1 Hz. Since earlier experiments had revealed 
troubles with clippers which did not give an equal number of zeros and ones, 

we introduced a 1-kHz offset in the local oscillator at Sugar Grove to avoid 

clipper bias at such a low fringe rate. This resulted in an artifical fringe 
rate of about 1.001 kHz on the short baseline and also added 1 kHz to the 
natural fringe rate on one long baseline (Arecibo to Sugar Grove) • Another 
problem with the short baseline was that the same interference signals might 
occur at both NRAO and Sugar Grove, and we could observe correlation due to 
interference. However, such a correlation would occur at the 1-kHz fringe 
rate and only the true fringes would show up at a fringe rate near 1.001 kHz. 

The Mark I recording system consists of a standard Ampex computer tape 

drive and a controller. The IF signal from the telescope is converted to 
video with a bandwidth of about 350 kHz. The controller clips the signal, 
samples it at 720 kbits/sec and sends the data, buffered in 6-bit characters, 
to the tape drive. One function of the controller is to start the tape re- 
cording at the time set in the schedule so that the run starts simultaneously 
with the other stations. The tape drive uses standard half-inch magnetic 
tape; it is run at its fast limit of 150 in/sec, 800 bits/ in. One record 
is 0.2 sec long, and there are about 900 records on one tape. Table A-2 
contains the data acquisition specifications. Because of the many reels 
of tape involved, the observing program at each station requires an obser- 
ver (tape-hanger) present at all times. More detailed information on the 
recording system is given by Hinteregger (1972). 
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Table A-2 

Mark 1 data acquisition specifications 


Quantity 

Bandwidth 

Sample rate = tape bit rate 

One bit 

One record 

Inter-record gap 

One tape 

One observation 

Tape speed 

Tape density 


Value 

350 kHz 
720 kbit/sec 
1.3888889 sec 
0.2 sec 
~5 msec 

900 records ■ 180 sec 

one tape 

150 Inches/sec 

800 6-bit characters/inch 
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In order to be able to synchronize the tapes during play-back In the 
computer we need to know the offset from UT of the VLBI master clock. We 
attempt to Insure that the clock at each station Is set to within a few 
microseconds of UT. This Is usually done In one of three ways: (1) We 

believe that the observatory clock has the correct time and use that. 

This has proven to be a disastrous assumption on enough occasions that 
(2) we require that the station clock be checked against Loran and WWV 
transmissions. For best Insurance, however, (3) we carry a rubidium clock 
standard to the Naval Observatory or the Goddard Space Flight Center clock 
laboratory, set the time, and transport the clock to the telescope. The 
clocks are "set" most accurately after the experiment by processing several 
tapes which were recorded on strong sources and searching for the clock off- 
set which yields the largest correlation coefficient. 




APPENDIX B 


INTERFEROMETER PHASE 


The expression used to calculate the a priori phase is derived below, 
and the result is used in the description of the operation of the two main 
processing programs, VLBI1 and VLBI2, in Chapter II, 

In general, the fringe phase is given by 

<Kt) * &>T g (t) + Aal Lo At LO (B -1 ) 

where u) = total local oscillator frequency, 

t = true geometrical delay which changes with hour angle, 

8 

= offset between local oscillators, 

t L Q = time elapsed since started at x^. 

Note that only the total local oscillator frequency enters the calculation. 

This is because the geometrical time delay, x , occurs in the RF signal and 

8 

the compensating delay introduced by the processing programs is inserted 
into the video signal. 

A Tay or series expansion is made for the true time delay: 

yo ■ yy + y t)(t - y (B - 2) 

where t = nominal start time of the observation. In universal time, 
o * 

t ■ time rate of change of geometrical delay. 

6 

The total calculated phase is 

* c .ic (t) ■ [“ r g ct o + y] + [<y°“ + 4 v “r] 

+ [ 4 “LO At Lo] <E_3) 


RECEDING PAGE BLANK 
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where x^ ■ difference between clock epochs r.t the two stations, 

Ax * time elapsed since the reference time t . 

K O 

Each term of equation B-3 will be discussed separately next. 

The first term in square brackets of eouation B-3 Is called the "funda- 
mertal phases" <|> o , 


I'T (t + T ) 

g o e 


(B-4 ) 


which accounts for the geometrical delays. The. delay is calculated at the 
nominal start of the observations, incl iding the known clock offsets. The 
next term refers to the change of j>.ionetrical lay as the earth rotates; 
thi3 is the fringe rate. The fringe phase is counted from the refer* nee 
time t Q . The "total fringe rate," f, is 


f - t w + Ao» . (B-5) 

g ID 

The last term is the phase of the oscillator which is used to offset the 
fringe rate if necessary. Its phase is counted from the time at which 
the oscillator was started. This time, t^, is specified in UT. The 
total calculated phase can now be written as: 

Ht) calc “ *o + fAt k + *“L0 At L0 ’ (B ~ 6) 

This expression is used in Chapter II, section B, as equation II-l. 



APPENDIX C 


REPLICATION OP PULSAR FRINGES 

Pulsar fringes are replicated because the source signal Is a periodic 
function. As an aid to understanding how this happens, figure C-l Is pre- 
sented. The symbols and functioi\al representations used here follow the 
style of Bracewell (1965). Explanatio s for the symbols and values of the 
numerical quantities used are given in table C-l. (These symbols are used 
only In this appendix and -re no' the same as in the main text.) 

Part (a) of figure C-l is a sample of a span of data recorded on a 
pulsar with s. period of about 0.62 sec using the Mark I system. An ana - 
lytic representation of the interferometer response as a function of time 
is shown in (b). Graphical sketches of each term are in (c). The fringe 
rate spectrum is the Fourier transform of the response; a representation 
is shown in (d). The convolution theorem for Fourier transforms was em- 
ployed in passing from the time domain to the frequency domain. Sketches 
ot each term in the fringe rate spectrum are given in (e) and the total 
spectrum is drawn in (f). The spectrum illustrated is a cut through the 
(tg, fg) plane of figure II-7 (page *.5) along the sloping axis. For most 
pulsars, the parameter G is very small so the spectrum lies nearly along 
the axis. 

The basic sampling interval, with the gating techrlque, is the pulbe 
period, and the noise spectrum extend' n‘.y out to frequencies ±1/2P Hr 
instead of t2. 5 Hz as for a continuum source. However. rote that the 
amplitudes of the Fourier harmonics are modulated Loth by the pulse sl.-pe 
transform and by the basic sampling window transform of full width i.0 Hz. 

The true amplitudes can be obtained by first demodulating with the calculated 
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REPLICATION OF PULSAR FRIN6ES 

— ji — ■ i i-j- L t .A. 


I L 


START i o 
TAPE LX- 
R 


-w 




TIME 


(a) SAMPLE SPAN OF DATA. 


R(t) = [Aft) * TT (T)] * [l P (t)*TT(W 6 )»lll(P)] * [lU(R)*W(D)]*[lll(R)] 


FRINGES OBSERVING 

WINDOW 


PULSE GATE PULSE RECORD DATA SAMPLING 

SHAPE WMOOW REPETITION REPETITION WINDOW FUNCTION 


(b) REPRESENTATION OF (a) IN FUNCTIONAL FORM. 


^ Tj *g P 

(c) GRAPHICAL REPRESENTATION OF (b) 


^ rilL-XLl.riII 




*- r R 


IP(f) = F(f)*sinc 0/t) * [F P (f)«*sinc 0/w 6 ) x m(i/p) * [ill (Vn)* sine 0/b)]« 1U( i /r)] 


FRINGE RESOLUTION 
SPECTRUM MODU- 
LATION 


PULSE GATE SAMPLING SAMPLING FRINGE REPLICATION 

SHAPE MODULATION FUNCTION FUNCTION MODU- OF PATTERN 

TRANSFORM LATION 


(d) FOURIER TRANSFORM OF (b). 
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1 rv-- , 

« 

jtllll k, 

* 

r 'T' 1 

/ 

/ \ 

i L 

M 

mi 
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(e) GRAPHICAL REPRESENTATION OF (d) 



(f) TOTAL FRINGE RATE SPECTRUM (Entire pattern replicated every Vb Hz). 


Figure C-l. Replication of pulsar fringes. Refer to table C-l 
and the text for explanation of the symbols used. The fringe 
amplitude in (f) will be modulated by the long-dashed curve as 
well as by the short-dashed curve. 
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Symbol 

t 

f 

P 


V 


p 


D 

I 


R 

T 

n(x) 

TTT(y) 


3inc(z) 

x 

* 


Table C-l 

Explanation of symbols in figure C-l 

Explanation 

time 

fringe frequency 

pulse period, M).62 sec in this example 
pulse width 

gate width, ~50 msec in this example 

data span ~190 msec 

inter-record gap ^5 msec 

record interval = 0.2 sec 

total length of observation = 180 ^ec 

box-car function of width x 

picket fence function: infinite series of delta functions 

spaced y apart 

sin(^z)/irz function of half-width z 

multiplication 

convolution 
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value of sine x at the appropriate frequencies. We could then turn to 
the Fourier component method described in Chapter II se. .ion B-3 and re- 
construct a pulse shape using equation 11-14. 



APPENDIX D 


INTERFEROMETER INSTRUMENTAL POLARIZATION 


A. Theory 

An unwanted response results from incomplete isolation of the right 
circularly polarized (RCP) feed from the left circularly polarized (LCP) 
feed of an antenna. The RCP feed responds to a fraction (f) of the LCP 
component of the signal, and only responds to a fraction (1-f) of the 
RCP component. The quantity f is the single- telescope instrumental polar- 
ization. The response of an interferometer formed by two telescopes which 
have non-zero values for f is reduced because the RCP and LCP components 
of the signal are not correlated. This effect can be corrected for if 
values of f, or some other calibrated quantity, are known. The expression 
for a correction factor is derived as follows. 

The total power observed with an RCP feed at stations 1 and 2 is. 


P 1R “ v f l )P lR + f l P lL 


2R 


(l- f 2> p 2 R + f 2 P 2L 


where and P££ * true RCP component of power, 

P°£ and P°^ * observed RCP component of power, 

f^ and f 2 “ fractional response of antenna to oppositely polarized 
radiation; "instrumental polarization" of single tele- 
scopes (ideally, fj * f^ “ 0). 

If we assume all the phase angle differences average to zero and there is 
no circular polarization in the source radiation, then ■ Pj^ and 

P 2R " P 2L* 
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Let 


e = £^2 - 2f x f 2 


and 


m “ P LR^ P RR 


The quantity m is the true fractional linear polarization of the source's 
correlated flux; and are the total and linearly polarized flux, 
respectively. When each antenna observes with an RCP feed the observed 
correlated power is 


P RR P RR 


1 - 6 + m 2 (3 


I 1/2 
6 


When one antenna observes with an RCP feed and the other with an LCP feed, 
the correlation coefficient is: 

1/2 


o I Q , 2 

P Rl. = P RR 1 6 + " 


6 + m 2 l 

- J 


Let the fraction a be the ratio of the correlation coefficient ob- 
served with similarly- and oppositely-polarized feeds and call it the 

"instrumental polarization" of the interferometer. Then a is a directly 

2 

observable quantity (we assume the term m g is negligible) : 

1/2 


o . o 

“ * V p i* 


>n 2 

L i - s 


The true correlation coefficient for a\ unpolarized source can be obtained 
bv: 


P RR - P RR (1 + “ > 


1/2 
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To determine m, an unpolarized calibrator source must be observed 

to obtain values of for the calibrator source. The polarized source 

determines a . Then 
P 


m = 


2 2 




1/2 


In principle, very low values of m can be measured even though is large, 

provided that and can be accurately determined. Figure D-l contains 

curves from which m can be read once a and a are known. 

P c 

Unless the source is circularly polarized we should observe that 

p = p TT and unless there is a nonlinear effect at one antenna we should 
HR LL 

also observe that P„ T * p 7D . Deviations from these equalities are attrib- 
utable to noise, coupling between dipoles, interaction between the feed 
and reed-support legs, or other unknown factors. 


B. Application 

The values of a were obtained from observations of strong extra- 
galactic sources which were assumed to be unpolarized. The sources we 
used were: 3C273, 3C286, 3C287, 3C298. At the NRAO station, the received 

signal was switched at a 2.5 Hz rate so that alternate records on the tape 
were recorded with the RCP and LCP feeds. At Arecibo and Sugar Grove, the 
entire tape was recorded with either the RCP or LCP feed. Then the cross- 
correlation output from the Arecibo-NRAO or NRAO-Sugar Grove baselines is 
alternately p_ D and p DT , or whatever combination occurred. An estimate 
of the quantity a can then be computed from each tape without the need 
for normalization and its attendant errors. The only assumption necessary 
is that the LCP and RCP feeds have the same efficiency. 
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This procedure was followed in 1972 January and February and the 
calibrations were assumed to apply for the remaining months since the 
feeds were the same. In 1972 December and 1973 February extensive polar- 
ization observations were made for calibration purposes. Since the feeds 
were different for each experiment the values of f and o are different 
also. The observed values of a and the derived values of f are in Tables 11-5 
and II-6 of Chapter II (page 39) . The calibrations and effect on the experi- 
ment are discussed in the main text (pages 37-40). 



APPENDIX E 


CATALOGUE OF FOURIER COMPONENTS AND PULSE SHAPES 
FOR THE CRAB NEBULA PULSAR AT 196 MHZ 

The tables which follow contain the amplitudes and phases of the 
complex correlation coefficients detected in each three-minute observa- 
tion of the Crab nebula pulsar at 196 MHz. The title of each table indi- 
cates whether the data is from short- or long-baseline runs. Each column 
of the tables is labelled as to its contents, a description of which follows. 

Columns 1 and 2 give the run number and date of each observation. 

Column 3 (TIME) contains the UT at which the observation started. Column 4 
(BASE) indicates the baseline on which the observation was made, and each 
letter indicates the station (A for Arecibo, N for NRAO, S for Sugar Grove)* 
Column 5 contains the TOTAL FLUX (in Janskys); this is the total time-averaged 
flux, i.e. the compact source flux. These numbers were used to compute the 
average fluxes presented in table III-l (page 47). Columns 5, 6, 8, 10, 12, 
and 14 are the amplitudes, A^ where n=0 to 5, and columns 7, 9, 11, 13, and 
15 are the phases, where n=l to 5, of the Fourier components of the aver- 
age pulse shape (see equation II-5, page 22). The top line for each run 
contains the components for n > 0, the second line is for n < 0. The phases 
are given in degrees relative to the phase of the n*0 component. The phase 
convention used for this table is that, apart from noise, phases of a pair 
of components for the same value of n should be equal. Runs for which only 
the amplitude of the n*0 component is listed were usually too weak, contained 
excessive tape errors after many processings, or had components which were 
too noisy to be reliably detected. 
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The fifth column from the right (3 SIG LEVEL) contains the 3a defec- 
tion level; any component amplitudes below this value were probably not 
detected. The third and fourth columns from the right (PHASE) indicate 
the signed deviation of the +n and -n component phases from equality and 
the rms deviation. These numbers were used in calculating the averages 
found in table IV-5 (page 81) . The second column from the right (PULSE PEAK) 
gives the height of the maximum of the pulse profile. This was used for 
editing purposes. The last column (FRACTION PULSING) indicates the fraction 
of the total flux which is contained in pulsing power; at 196 MHz all of 
the power is pulsing. 

Plots of the pulse profiles follow the tables. The baseline, fre- 
quency, and the month and year are indicated; the flux scales are in 
Janskys. The profiles have been superimposed by simply matching the peaks. 
The average profiles are therefore incoherent averages. 
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Figure E~l. Individual and average pulse shapes at 196 MHz, 





Figure E-l. (continued) 



APPENDIX F 


CATALOGUE OF FOURIER COMPONENTS AND PULSE SHAPES 
ioR THE CRAB NEBULA PULSAR AT 111 MHZ 


The tables which follow contain the amplitudes and phases of the 
complex correlation coefficients detected in each three-minute observa- 
tion of the Crab nebula pulsar at 111 MHz, Explanations of the contents 
of each column were given on the f irst page of Appendix E, page 150. 

For the 111-MHz data, the last column (FRACTION PULSING) was used to 
compute the averages found in table IV-3 (page 78) and plotted in fig- 
ure IV- 5 (page 77). Although three Fourier components are listed for 
each run, the n = ±3 components are often not detected in December 1972 
and February 1973. For these dates, only the first two components were 
used in constructing the pulses. 
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